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Abstract 

We introduce an individual-based model for fiber elements having the ability 
to cross-link or unlink each other and to align with each other at the cross links. 

We first formally derive a kinetic model for the fiber and cross-links distribution 
functions. We then consider the fast linking/unlinking regime in which the model 
can be reduced to the fiber distribution function only and investigate its diffusion 
limit. The resulting macroscopic model consists of a system of nonlinear diffusion 
equations for the fiber density and mean orientation. In the case of a homogeneous 
fiber density, we show that the model is elliptic. 
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1 Introduction 

The topic of complex systems is attracting an increasingly abundant literature, due to 
its paramount importance in life and social sciences. Complex systems consist of a large 
number of agents interacting through local interactions only and yet able to self-organize 
into large-scale coherent structures and collective motion [36]. Among examples of in¬ 
teractions leading to collective motion, the alignment interaction has been the subject of 
many studies since the seminal work of Vicsek and co-authors [35]. In Vicsek’s model, 
self-propelled point particles tend to align with their neighbors up to some noise. Vicsek’s 
particles are polar: they carry a definite direction and orientation defined by the unit 
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vector of their propulsion velocity. Their alignment interaction is also polar in the sense 
that a particle moving in an opposite direction to its neighbors will eventually reverse 
its direction of motion. However, other alignment rules have been studied as well. Polar 
particles can be subjected to nematic alignment. In this case, a particle moving in an 
opposite direction to its neighbors will not reverse its direction of motion, as opposed to 
the polar alignment case. Nematic alignment has been used as a model for the volume 
exclusion interaction [5j [20(1291 .Particles can also be apolar, for instance if they randomly 
reverse their direction of motion. Apolar particles interacting through nematic alignment 
have been proposed as a model for vibrating rods [6], or fiber networks pQ . In the related 
held of nematic liquid crystals, volume exclusion interactions between rod-like particles 
are also modelled as an alignment force [TBj [24], [28]. But additionally, the molecules are 
convected by the background solvent and are subjected to rotation by the fluid shear. 
Additionally, they contribute to the fluid dynamics of the liquid solvent through an ad¬ 
ditional extra-stress tensor. Usually, the polymer chains are supposed of fixed length, 
although lately, models of polymer chains of variables lengths have appeared [T2] . 

In the present work, we are interested in a system consisting of fibers (or polymer 
chains) of variable lengths. This model aims to describe the network of collagen fibers in 
a fibrous tissue. We model fiber length variation (through polymerization / depolymer¬ 
ization) as well as the ability for the fibers to establish cross-links between them by the 
same basic rules described as follows. We assume the existence of a fiber unit element 
(or monomer) modeled as a line segment of fixed length L. We suppose that two fiber 
elements that cross each-other may form a link, thereby creating a longer fiber. There 
is no limit to the number of cross-links a given fiber can make. Therefore, the fibers 
have the ability to branch off and to achieve complex network topologies. We include 
fiber resistance to bending by assuming the existence of torque which, in the absence of 
any other force, makes the two linked fiber elements align with each other. Fibers are 
also subject to random positional and orientational noise and to external positional and 
orientational potential forces. Finally, cross-links may also be removed to model possible 
fiber breakage or depolymerization. 

Our model features apolar fiber particles (since they are not self-propelled), interacting 
through nematic alignment with the other fibers they are linked to. Thus, the model bears 
analogies with previous models of apolar particles interacting through nematic alignment 
M. However, the interaction network topology (which keeps track of which fiber pairs 
are cross-linked) is different, as ours is determined by the distribution of cross-links. 
The fact that this network topology changes with time through dynamic cross-linking or 
unlinking processes is one specific feature of the present work. In the absence of cross¬ 
link remodeling, i.e. when the cross-links lifetime is infinite and no new cross-links is 
created, each connected component of the fiber network can be seen as an unstretchable 
elastic string since all connected fiber elements will spontaneously align with each other. 
However, cross-link removal or creation events (supposed to occur at Poisson distributed 
random times) introduce a fluid-like component to the rheology of the fibers, thereby 
confering some visco-elastic character to the medium. Cross-link-governed statics and 
dynamics of fiber networks have been intensely studied in the literature [3j [8, :9| (21, [27] 

. However, most models consider passive cross-links which only act on the fibers by a 
spring-like attractive force. Here, our description introduces active links which tend to 


2 


align the two fibers with each other. By doing so, we are also able to take into account 
fiber breakage, elongation and branching just in addition to and in the same way as fiber 
linking/unlinking because cross-linked fiber elements can be seen as two parts of the same 
fiber. Another difference from previous literature is that fibers in our model are subject 
to noise making the system more akin to a fluid or a gas than to a solid. By contrast to 
classical polymeric fluid studies, we do not assume that the fibers are transported by a 
fluid and modify its rheological properties but this feature could be added in future work. 

This model was first introduced in Ref. [30] where it was coupled with the dynamics 
of spherical particles modelling cells. This model has been built to describe the self¬ 
organization of the adipose tissue, where spheres represent adipocytes and fibers, the 
surrounding collagen fibers. In this work, we demonstrated that the interaction between 
cells and fibers led to the spontaneous formation of cell clusters of ovoid shape akin to 
the adipose lobules that form the functional subunits of the adipose tissue. In Ref. [30], 
only a discrete Individual-Based Model (IBM) was considered. The present work focuses 
on the fibrous medium only and aims to derive meso and macroscopic models from the 
background IBM using techniques of kinetic theory. Indeed, the computational cost of 
an IBM scales polynomially with the number of agents, which makes them practically 
untractable for large systems. Continuum models allow to break this curse of scaling but 
they suppose that a suitable coarse-graining procedure which averages out the fine-scale 
structure has been applied to the IBM. In order to capture the correct effects of the fine- 
scale dynamics on the large-scale structures, it is of paramount importance to perform 
this coarse-graining as rigorously as possible. This is the aim of the present work. 

The derivation of a continuum model from the fiber dynamics is done in two steps. 
We first derive a kinetic model from the underlying IBM and secondly, we perform a 
diffusion approximation of the latter to obtain the continuum model. The kinetic model 
provides a statistical mechanics description of the underlying IBM by investigating how 
the probability distribution of fibers in position and orientation space evolves in time. 
Here, we will show that the mere distribution of fibers is not sufficient to close the system 
and that the cross-link probability distribution needs to be introduced. The cross-links 
provide correlations between the fibers and consequently their distribution can be viewed 
as similar to the two-particle fiber distribution. We will formally show that the knowledge 
of the one- and two-particle distributions is enough to provide a valid kinetic description 
of the system. Of course, this fact needs to be confirmed by numerical simulations and 
mathematical proofs. But if it proves correct, this model provides a unique example, to our 
knowledge, of a kinetic model which is closed at the level of the two-particle distribution 
function. Indeed, the question whether or not kinetic descriptions must include higher 
order distribution functions has been actively discussed in the recent years [mi im 1251126] 

. We also note that the introduction of the cross-link distribution functions provides 
an economic and efficient way of statistically tracking the fiber network topology. This 
methodology could prove interesting for other situations of dynamically evolving networks. 

The second step consists of a diffusion approximation of the previously derived kinetic 
model. It starts with changing the time and space units to macroscopic ones. The 
macroscopic space unit is large compared to the typical spatial scale of the fibers, e.g. their 
length and the macroscopic time unit is large to the typical time scale of the fibers, e.g. 
the time needed for two linked fibers to align with each other. A diffusive rescaling relates 
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the time and space rescaling in such a way that the ratio of the microscopic to macroscopic 
time units is the square of that of the spatial units. This choice is made necessary by the 
absence of any polarization in the medium which makes diffusive behavior dominate. A 
key assumption that we make here is to assume that the linking/unlinking frequencies are 
very large: the typical linking/unlinking time measured in the macroscopic time unit scales 
like the square of the typical fiber alignment time (also measured in macroscopic unit), 
which is very small. This allows us to deduce an algebraic relation between the cross-link 
distribution function and the fiber distribution function, and to realize a closure of the 
kinetic equation at the level of the fiber distribution function alone. This assumption is 
questionable given the biological applications we have in mind, but it provides a first step 
towards a more complete theory involving finite linking/unlinking times. 

From these assumptions, we derive a singular perturbation problem for the fiber kinetic 
distribution function that has the form of a classical diffusion approximation problem 
® dm nn], whose leading order collision operator comes from the nematic alignment of 
the fibers due to the alignment torque at the cross-links. This operator has equilibria 
in the form of generalized von Mises distributions of the fiber directions. The von Mises 
distribution extends Gaussian distributions to probabilities defined on the unit circle. It 
is peaked around a mean fiber direction angle 6 0 . The continuum model describes how the 
local fiber density p and the local fiber direction 6 0 vary as functions of position x and time 
t. To obtain these evolution equations, we must integrate the kinetic equation against 
suitably chosen collision invariants. This operation cancels the singularly perturbed term. 
Here, the difficulty it that there exists only one such collision invariant in the classical 
sense, which allows us to find an equation for the density p only. To find an equation for 
the mean fiber direction 9 0l we use the recently developed theory of Generalized Collision 
Invariants (GCI) [TT1 [151 [171 IT9] . The resulting system is a nonlinear coupled system of 
diffusion equations for p and 6q. In the case of a homogeneous fiber distribution, when 
the density is uniform in space and constant in time, we show that the resulting nonlinear 
diffusion model for 6 0 is parabolic. In future work, it will be shown that this system is 
well-posed. Numerical simulations will demonstrate that the continuum model provides 
a consistent approximation of the underlying IBM for the fiber dynamics. Numerous 
macroscopic models for fibrous media have been previously considered in the literature 
but very few of them have been derived from an underlying IBM. Most of them are 
heuristically derived from continuum theories such as mechano-chemical principles [2} S3] , 
thermodynamics [22], or viscous fluid mechanics [23]. 

The outline of this paper is as follows. In Section [21 we start with the description of the 
IBM. Section [3] is devoted to the derivation of the kinetic model. The scaling assumptions 
and the scaled kinetic equations are derived in Section [4j In Section [5j we perform the 
large scale limit of the so-obtained equations. Finally, Section [6] is devoted to the analysis 
of the model in the case of a homogeneous fiber density. Conclusions and perspectives 
are drawn in Section [71 Some technical computations are detailed in Appendices. 
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2 Individual Based Model for fibers interacting through 
alignment interactions 

We intend to model a medium consisting of interconnected fibers. To simplify the ge¬ 
ometric description of fibers, we decompose them into fiber elements of uniform fixed 
length and consider that a fiber consists of several connected fiber elements. The link 
between two connected fibers can be positionned at any point along the fibers (not only 
the extremities) and a given fiber can be connected to any number of other fibers, thereby 
allowing to model the branching off of a fiber into several branches. The links are not 
permanent. The topology of the fiber network is constantly remodelled through link cre¬ 
ation/deletion processes. To model fiber resistance to bending, we suppose that pairs 
of linked fibrs are subject to a torque that tends to align the two fibers with respect to 
each other. Finally, the fibers are subject to random positional and orientational noises 
to model the movements of the tissue and to positional and orientational potential forces 
to model the action of external elements. In the case of a fibrous tissue, these external 
elements may consist of cells or other tissues. 

In this paper, we restrict ourselves to a two-dimensional model. We consider a set of 
N fiber elements modelled as small line segments of uniform and fixed length L , described 
by their center X 7 6 M 2 and their angle 0 7 with respect to a fixed reference direction. As 
the fiber elements are assumed apolar, 0 7 is an angle of lines, i.e. 9 7 G [—|, modulo ^ r. 

We define energies related to each of the phenomena described above namely an energy 
for the maintenance of the links Winks, an energy for the alignment torque IF a i ign , an 
energy for the action of the external elements W ext , an energy for the noise contribution 
Woise and a total energy made of the sum of all these energies: 

Wtot = Winks + W ex t + IFalign + Wioise, (2-1) 

All these energies are functions of the N fiber positions (A" i )^ 1 and orientations (9i)^ =1 . 
Note that fF n oise is rather an entropy than an energy, so that W tot is indeed the total free 
energy of the system. Fiber motion and rotation during a time interval between two fiber 
linking-unlinking events is supposed to occur in the steepest descent direction to this free 
energy, namely according to: 


dXi _ 
dt 

—h VjsqlFtot, 



(2.2) 

(19 7 
dt 

—A d g W tot , 

Vi e {!,.. 


(2.3) 


Eqs. (12.211 and (12.31) express the motion and rotation of the individuals in an overdamped 
regime in which the forces due to friction are very large compared to the inertial forces. 
Fiber velocity and angular speed are proportional to the force exerted on the fiber through 
two mobility coefficients /i and A which are considered given. We now detail the expres¬ 
sions of the four energies involved in the expression (12.11) of the total free energy of the 
system, as well as how Eqs. (12.21) and (12.31) are supplemented by Poisson jump processes 
when a linking/unlinking event occurs. 

To define the expression of Winks, we consider a time at which no linking/unlinking 
process occurs. Then, the set of links is well-defined and supposed to have K elements. Let 
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k G {1,..., K} be a given link and denote by ( i(k),j(k )) the pair of indices corresponding 
to the two fibers connected by this link. To make the labeling of the pair unique, we assume 
without loss of generality that the first element of the linked pair is always the one with 
lowest index, i.e. i{k) < j{k). The link is supposed to connect two points XW and 
X*X on fibers i{k ) and j{k) respectively. These points are determined by the algebraic 
distances £ k ^ and to the centers X^ and Xj^ of the two fibers respectively; We 
thus have the relation: 


X i{k) -Xi(k ) + X j(k) 

where £ k t k y G [— L/2,L/2] and where, for any fiber i, we let a \ = (cos 9i, sin Oj) be 
the unit vector in the direction of the fiber. All along the link lifetime, the link places 
a spring-like restoring force that attracts X^ back to Xj^) (and vice-versa) as soon as 
their are displaced one with respect to each other. This restoring force gives rise to a 
potential energy R(X i(fc) , £^ k) , X j(fc) , 6> j(fc) , £ k j(k) ), with 

V (X 1? 9i,£i,X 2 ,9 2 ,£ 2 ) = -\X\-\- £\uj{9i) — (X 2 + £ 2 uj(9 2 ))\ 2 , (2-4) 

where k, is the intensity of the restoring force. Obviously, the larger k, the better the 
maintainance of the link is ensured. The potential Wii n > s is then assumed to be the sum 
of all the linked fiber spring forces: 


R links 


1 

2 


K 

Y. X m , 

k =1 


(2.5) 


We stress the fact that the quantities £ k ^ and £ k -^ remain constant throughout the link 
lifetime. They are determined at the time of the creation of the link (see below and 

Fig. CD). 



Figure 1: Intersecting linked fibers. kj and Iji refer to £(Xi,9i, Xj,9j) and 
£(Xj, 9j, Xi, 9i) (12.121) . A. Situation at linking time. B. Restoring potential V l3 (12.41) 
after motion of the fibers. 


The external potential W ext associated with the external forces is supposed to be the 
sum of potential forces f/(Xj, 9i ) acting on each of the N fibers: 

N 

W ext = J2 u ( x ii9i). ( 2 . 6 ) 

i=l 


6 







Here, U(x,9) is a given, possibly time-dependent smooth function. In the case where the 
system describes the collagen fibers in a tissue, U aims to model the presence of cells or 
other organs. 

Linked fibers are subjected to an alignment force at their junction to model fiber 
resistance to bending. This force tends to align linked fibers i(k) and j(k) and derives 
from the potential &( 6 ^(fcy 9j^) which reads: 

b{e u e 2 ) = a\ sin(0 1 -0 2 )|' J , (2.7) 

where a plays the role of a flexural modulus and /3 is a modeling parameter. The binary 
alignment potential only depends on the angles 9± and 9 2 , and the total alignment energy 
Waiign is supposed to be the sum of all the binary alignment interactions: 

1 K 

Waiign = b(9i( k ), 9j(k)). (2.8) 

z k =1 

We include random positional and orientational motion of the fiber elements which, 
in the context of tissue dynamics, originate from the random movements of the subject. 
With this aim, we introduce an entropy term: 

N 

W noise = d^og(f)(Xi,9i), (2.9) 

i =1 


where / is a ’regularized density’: 


1 N 

f{x,0) = jj'52£ N ( x - x i)v N (0-0i)- 

i =1 

Here, £ N and r] N are regularization functions which allow to define the logarithm of / and 
have the following properties: 

'fe^(l 2 ), V N £C~([-^}), £ n > 0, r] N > 0, 

' J(. N (x)dx = 1 , J = 1 , 

.Supp(y) C B(0,R"), Supp(^) c [ -M n ,M n ], 

where C'°°(M 2 ) is the set of infinitely differentiable functions on M 2 , C^ r ([—|, |]) the set 
of periodic C°° functions of [—|, |] and Supp stands for the support of a function. Here, 
R N and M N are chosen such that and NM n — y oo as —y oo. The mean 

interparticle distance in x and 6 are respectively of order - 7 = and -L. This condition is 
equivalent to — > 0 and N ^ fN —> 0, which means that as iV —> 00 , the number of 

particles inside the support of a regularizing kernel tends to infinity. This way of modeling 
the influence of the noise is customary in polymer dynamics jT] .In the next section, we 
show that such an entropy term gives rise to diffusion terms at the level of the mean-held 
kinetic model. 
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By inserting (j2.5|) . (12.61) . (12. 8 p and (12.9p into (|2.2j) . (12.3p . we find the fiber equation of 
motion, during any time interval between two linking/unlinking events: 


dXj 

dt 


= - V 


w x (u + \ogf N )(x. l ,e l ) 


K 


d" r) "y 1 V Xl V (Xj(fc), Xj(fc), £/( fc )) 

k=l,i(k)=i 
K 

y i ^ x 2 vpi(k)i q j(k)i ^j{k )) 

k=l,j(k)=i 


2 

1 


d(h 

dt 


A 

1 

2 

1 

2 

1 

2 


+ log/ JV )(X i , 0j) 

K 

£ a»,v(x i(k) ,e m ,e k lk) ,x m ,e m ,e k (k) ) 

k=l,i(k)=i 

K 

£ A >>. 

fc=l 5< 7(/c)=i 

K 1 K 

y 1 ^j(fc)) + 2 /* 1 @q 2 b{6i(k),6j(k)) 

k=l,i(k)=i k=l,j(k)=i 




which we can write: 


dt 


= -/i 


iC 


iC 


x V i i(t) (i)VJt-^ S m (i)X„V (C\ 


k=1 


fc=l 




+ V K (f/ + log f N )(Xi, 6i) 


ddi 

dt 


= - A 


w+iognpM,) 

X 1 K 

2 ^(fc)(*)^ 0 i^ + 2 ^'( fc )(*)^ e 2 ^ ) {Ci(h),j(k)) 

k= 1 /c=l 

1 X x X 

2 1 ^ ^ 8j(k)(i)dg 2 b J (dj(/e), dj(k)) 

k= 1 /c=l 


( 2 . 10 ) 


( 2 , 11 ) 


with C k k)J(t) = (X i(k) ,ff i ( k) ,e k (k) ,X j f k ),6j( k ),l k ( k) ) and <5 i(i!) (i) is the Kronecker symbol, i.e. 

= 1 if i(k) = i and S^k)(i) = 0 otherwise. 

When two fibers i and j intersect each other, because of the continuity of their motion, 
they are going to intersect each other during a time interval [t*, t*]. We assume that, during 
this time span, the linking probability follows a Poisson process of parameter i/f, i.e. the 
probability that a link is formed during the interval /*, t] with t < t* is 1 — 






















Only one link can be formed between the two fibers of the same fiber pair. Supposing 
that a link, indexed by k is formed between the fibers i and j (such that i — i(k) and 
3=m if i < j ) at a time t k G [£*,i*], we denote by X k the attachment site of the link. 
The distance £(Xi^,9i^,Xj^k),dj(k)) between the center X^ k ) of fiber i[k) to the k-th 
link attachment site X k with fiber j(k ) (see Figure HJB) can be directly computed by: 




rj'W 




sin(6> j(fc ) - 6>j ( k )) 




where = (x^, y^) are the coordinates of the center of fiber i(k). For X - 

and u = (ce,/3), we denote by X x u = x/3 - ya. Then, £(X^ k) , 

written: 


@i(k)i Xj{k)i ^j (k )) 


-W(fc)) X 0j(9j( k ')\ 


= (x,y) 

can be 


where again, uj(9) = (cos 9, sin 9) is the directional vector associated to angle 9. The fact 
that the two fibers are intersecting each other at time t k is written: 


K(-^Q(fc)j @i(k)i Xj(k)i @j(k)) I — 2’ all d 11 ®j{k) j 5 @i(k )) | ^ 

where L is the fiber length and where all positions and angles are evaluated at time t k . 
The quantities l(X i{k ), 9^ k ), ^j(fc), Oj{k)) and I(Xj( k ), 9j( k ), -W(fc), 9^ k )) at the time t k of the 
formation of the link set the positions of the attachment sites X^ and X k j(k) of the link 
on fibers i and j. Therefore, and £ k ^ remain constant throughout the link lifetime 
and equal to their value at the time t k . So, we have 


^ pk _ ^ pk _ n 

~ ~dt m ~ u ’ 


dt 


throughout the lifetime of the link. 

We also assume that existing links can disappear according to a Poisson random 
process of parameter i.e. the probability that the link disappears in the time interval 
[t k ,t] is l-e- Ud ^- tk \ 

The next section is devoted to the asymptotic limit N, K — y oo of this model. 


3 Derivation of a kinetic model 

Here, the derivation of a kinetic model from the Individual Based Model of section [2] is 
performed. The empirical measure f N (x,9,t) of the fibers is introduced: 

1 N 

f N (x, 9, t) = 

i= 1 

where S(x,xt),o,Xt))(x, 9) denotes the Dirac delta located at (Xj(t), 9j(t). It gives the proba¬ 
bility to find a fiber at point x and orientational angle 9 at time t. The empirical measure 
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g K (xi, 9\, £\, x 2 ,9 2 , £ 2 ,t) of the fiber links is given by: 


K 


g K (x 1 ,9 1 ,£ 1 ,x 2 ,9 2 ,£ 2l t) = 2 ^ aT2 > 


k=1 


+ S (x m ,e m /* (k) ,x iW A w ^ k) )( x ^ 0i,ti,x 2 ,9 2 ,e 2 ), 


with a similar definition of the Dirac deltas. It gives the probability of finding a link 
with associated lengths within a volume d£\d£ 2 about l\ and £ 2 , this link connecting a 
fiber located within a volume dx±^- about {x\,9\) with a fiber located within a volume 
dx 2 — about ( x 2 ,9 2 ). One notes that ( 7 , 7 ) is defined in [— k, ^] 2 . Then, at the limit 
N, K —> 00 , jf ~>■ 7 where £ > 0 is a fixed parameter, f N —$■ /, g K —> g where / and g 
satisfy equations given in the following theorem: 

Theorem 3.1. The formal limit of Eqs. (I2.2[ ) . (I2.3j) for K , iV — y 00 , —y where f > 0 

is a fixed parameter reads: 


df_ 

dt 

and 


-li[V x -((V x U)f)+fV x -F 1 + dA x f)-X[d e ((d e U)f)+fdgF 2 + dd 2 d f) =0, (3.1) 


dg 

dt 


-1* fa X1 ■ (gVMxM + Z ° F^xM) 

+ V X2 • (, gW x U(x 2 , e 2 ) + 9 2 )) 


9 


dy xi 'ir/ n 

f{x i,9i 


-V77iA)) + dV 


X2 


7 A, 9 2 


;Vxf(x 2 ,9 2 )) 


~x(d ei (gd e U(x l ,9 1 ) +£ ^ ^ F 2 (xi,6>i)) 

+ <% 2 (gdeU(x2,9 2 ) +Z X , 9 a , F 2 (x 2 ,9 2 )) 

J\ x 2 , v 2 ) 


dd 01 ( , 9 ^ d e f(x i,fli)) + dde 2 { £ ^ 9 , d e f(x 2j 9 2 )) ) = S^), 


77i7i) 


772 , 6 * 2 ) 


where 


d9 2 


Fi(x 1 , 6 * 1 ) = J (g\ r Xl V)(xi,9 1 ,£i,x 2 ,9 2 ,£ 2 )d£ 1 d£ 2 —dx 2 , 

/ d0 

(g(d 01 V + <%7)) 7i, 6 >i, 0, £ 2 , 6 * 2 , £ 2 )d£ 1 d£ 2 -^-dx 2 , 

and S(g) is given by: 

S(g ) = l'//(^lVl)/(^2,#2)% I1 ,8 1 , ra ,fe)(<l)y ra ,fe, I1 ,» 1 )(fe) - 


(3.2) 


(3.3) 

(3.4) 

(3.5) 


where 7(7) denotes the Dirac delta at £, i.e. the distribution acting on test functions 
0 (-7) such that (7(7), <77)) = <77 
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This kinetic model consists of two evolution equations. The first one (Eq. (13. ip ) is an 
equation for the individual fibers and describes the evolution of the one-particle distribu¬ 
tion function /. Eq. (I3.2j) is an equation for the links between fiber pairs. The distribution 
function g describes the correlations between fiber pairs brought by the presence of links. 
It can be viewed as a kind of two-particle fiber distribution function. This model is, to 
our knowledge, a unique explicit example of a kinetic model written in terms of the one 
and two particle ditributions and closed at this level. Also, the distribution function g 
can be seen as a way of describing the random graph of the fiber links, namely the graph 
where the nodes are the fibers and the edges are the links. This statistical description of 
a random graph could be useful to describe other kinds of random networks, notably in 
social sciences. As the links are tightly tied to the fibers, they are convected by them and 
follow their motion. Simultaneously, they constrain the linked fibers to move together, 
so they directly influence their motion. The action of the links on the individual fiber 
motion is contained in the third and sixth force terms F\ and F 2 of Eq. (13.ip and are the 
kinetic counterparts of (12.4p . The second and fith terms describe transport in physical 
and orientational spaces due to the external potential and are the kinetic counterparts 
of (12.6p . The fourth and seventh terms are diffusion terms of amplitude A d and /id respec¬ 
tively. They represent the random motion of the fibers and originate from the interactions 
described by Eq. (12.9p . The individual motion of the fibers is thus related to the motion 
of its linked neighbors. The left-hand side of Equation (13.2p describes the evolution of 
the links between fibers. Indeed, it is composed of the convective terms generated by the 
external potential and by the diffusion terms. The forces induced by the restoring poten¬ 
tial generated by the links again gives rise to the non local terms iq and the first term of 
F 2 . The kinetic counterpart of the alignment force between linked fibers (see Eq. (I2.8p l 
is encompassed in the second term of the force F 2 and only acts on the orientation of 
the fibers. The right hand side S(g) of equation (13. 2 p describes the Poisson processes 
of linking/unlinking at frequencies Uf and rq, respectively. The first term describes the 
formation of the link and the Dirac deltas indicate that, at the link creation time, the link 
lengths i\ and i 2 are set by the geometric configuration of the fibers at the attachment 
time. Also, because £\ and £ 2 are restricted to lie in the interval [— L/2,L/2], we see 
that the link creation term is non-zero only when two fiber elements are intersecting each 
other. The second term just describes a decay of the link distribution at the rate set by 
the Poisson process, i.e. u d . 

The formal proof of this result is inspired from Ref. |32j, and the detailed computations 
can be found in appendix [Aj The rigorous proof of this result is an open question and is 
left for future work. 


4 Scaling 

4.1 Dimensionless Equations 

We express the problem in dimensionless variables. For this purpose, let f 0 be the unit of 

1 1 X 2 

time and £o, /o = 73 -, 9o = rs and Uq = -A the units of space, distribution function and 

x 0 x 0 *0 

energy. The scaling of f(x, 6) and g(x\, 9i,£\, x 2 , 0 2 , £ 2 ) comes from the fact that they are 


11 













probability distribution functions on a 2D domain. The following dimensionless variables 
are defined: 


x 


x = 


f _f _2 - _ 9 _ fi T~T _ { 0 U 


= —, f = T = f x o> 9 = — = gx0, U = 


Xn 


xo X 0 fo ~ 9o 

and the following dimensionless parameters are introduced: 

li = —, A = y—, = t 0 u f , u d = t 0 v dl L = —, d = —a = —, /c = Kto- 


tn 


tn 


X 0 






First of all, from the expression of V (see Eq. (12.4ft h we get: 


V(x\, 9\, £ 1 , X 2 , 0 2 , £ 2 ) — + ^i^{9\) — x 2 — £2^{9 2 )) 2 

Zt 0 

— -jjV(Xl,Oi,£i,X2,02,^2), 

t 0 

with 

E(xi, 6 »i,li,x 2 , 6 l 2 ,l 2 ) = +^io;(6 l i) - x 2 - ^(6* 2 )) 2 . 

Now, from Eqs. fl3.3p ~ fl3.4p . one notes that: 


Fi(xi,6»i) = -^(xiA), 

where 

Ei(xi,6>i) = f V Sil V{x 1 ,9i,£i,x 2 ,9 2 ,£ 2 )g(xi,9 1 ,£i,x 2 ,9 2 ,£ 2 )d£id£ 2 —dx 2 , 

Jo 77 

with C! = R 2 x [—|, |] x [—y, y] x [—y, y]. Similarly, F 2 (a; 1 ,0 1 ) = ^F 2 (^i,^i), where: 
E 2 (xi, 6»i) = F aJ (xi, 6»i) + E /irifc (xi, 61 ), 

f - - d9 

Fnnk(x iA) = / (gde 1 V)(xi,9 1 ,£ l ,x 2 ,9 2 ,£ 2 )d£ l d£ 2 —dx 2 , 

Jo 77 

r _ _ (jCj 

F a i(x 1 ,9 1 )= / (gde 1 b)(x 1 ,9 1 J 1 ,x 2 ,9 2 J 2 )d£ 1 d£ 2 —dx 2 , 
k Jo 71 

where b(9i,9 2 ) = a/sin(0! — 02) /3 - In this new set of variables, Eqs. fl3.ip ~ fl3.2p become: ; 


dff-x W £ • (V*tf/) - X'deidgUf ) - eA'ckA-xCA'V* • A 

— d'X'dgf — d'xX'A x f = 0, 
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and 


dt'9-X A'Vsj • {gV s U(x i, 6> x ) + £ 7 ^ Fi(xi,0i)) 

f(xi,0 1) 

— A'd^ (gdeU(xi, 6>i) + £ r,® a 

fix i,9i) 

-XA'V, 2 • (flV,£/(x 2 ,0 2 ) + £ .J* . . FifeA)) 

f{x 2l 9 2 ) 

—\'do 2 (gdoU(x 2l 0 2 ) + £y—^——F 2 (x 2 , d 2 )) 

f{x 2 ,e 2 ) 

-d'x A'Vft • ( ,/, Vj(xM) 
fix i,9i) 

-d'x W a • ( .-, g V,/(i 2 ,g 2 )) 

fiX2,0 2 ) 

-d!\'d 0l i j.J* M {x i,fli)) 
fix i,9i) 

-d'\'d,,( j-?-— dJ(x 2 ,e 2 )) = S(g), 

J j y 2j 

where X = f 7 anc ^ : 

S(s)(&i, 0i, !i, x 2 ,0 2 , 4) =v' f Jix 1 ,0i)/(z 2 ,02)% 1 ,0 1 ,s 2 ,0 2 )(^i)%2,e2,si,ei)(4) 

- v' d g{x 1 ,d 1} £ 1 ,x 2 ,9 2 ,I 2 ). 

Finally, if the space and time scales xo, to are chosen such that A' = y = 1, i.e: 

Xq = —, to — Ah 

the dimensionless equations for / and g read (dropping the primes and tildes for the sake 
of clarity): 

'd t f - V x • (V,tt/) - deideUf) - id e F 2 - £V X • F, - d«9 2 / - dA x / = 0, (4.1) 

% - V X1 • (^V ai £/(x 1 ,0i) + £ 9 a . Fi(xi,0 X )) 

fix i,0i) 

-a ei (^C/(xi,0i) + £ ^ F 2 (xi,0i)) 

j(z 1,01) 

-V X2 • (gV x t/(x 2 ,0 2 ) + £ 9 a . Ffa, 0 2 )) 

f i x 2i 9 2 ) 

-de 2 igdeUix 2l 9 2 )+i 9 F 2 (x 2 ,0 2 )) (4.2) 

< jiX2,9 2 ) 

' { nfbrF J{xM) 

~ dd 0ii ff 9 a s d efix l,0l)) 

/(Zl,01) 

-dd^J^defM)) = S(g)(xi,9 1 ,£i,x 2 ,9 2 ,£ 2 ), 
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with 





J c 


F 2 (x 1: 9 1 ) = F al (x i, 6»i) + F link (x i, 0x), 



£ 





J c 


where C = M 2 x [— |, |] x [— L/2, L/2] x [—L/2, L/2] and 

ffi, £i, x 2 ,0 2 , i 2 ) =ie f f(x 1 ,9 1 )f(x 2 , 0 2 )Se( Xl ,e 1 ,x 2 ,e 2 ){ii)^(x 2 ,9 2 ,x 1 ,e 1 )(h) 

- u d g(x 1 ,9 1 J 1 ,x 2 ,9 2 ,£ 2 ). 

4.2 Scaled equations 

So far, the chosen time and space scales are microscopic ones, and describe the system 
at the scale of the agent interactions. In order to describe the system at a macroscopic 
scale, a small parameter e < 1 is introduced and the space and time units are set to 
£q = £ -1 / 2 £o, to — £ -1 h)- The fiber length measured at scale x 0 is supposed to stay 
of order 1 as e —> 0, i.e. L = 0(1). The variables x, t, i and unknowns / and g are 
then correspondingly changed to x — yj ~ex , t — et, £ — sjel, f(x,9 ) = e~ 1 f(x,8 ) and 
g(x\, 9i,£i, x 2 , 9 2 , £ 2 ) = e~ 3 g(xi, 9i, £±, x 2 , 9 2 , i 2 ). We suppose that the external potential 
U(x,9 ) is decomposed into U(x,9) = U°(x) + 0 1 (6 I ), where U° is acting on the space 
variable only and U 1 is a 7r-periodic potential acting on fiber orientation angles only. The 
external potential acting on the space variables is supposed to be one order of magnitude 
stronger than the one acting on the fiber rotations: U° = 0(1), U 1 = 0(e), i.e. U 1 = 
e^U 1 with U 1 = 0(1). The strength of the alignment potential is supposed to be large 
a = 0(s" 1 ), i.e. a = ea with a = 0(1), and we choose the exposant /3 = 1. The intensity 
of the alignment potential between linked fibers is supposed to be small k, = 0(e), i.e. 
h = e _1 K with k = 0(1) and the diffusion coefficient and parameter £ are supposed to 
stay of order 1: d,£ = 0(1). In order to simplify the analysis of the system, the process 
of linking/unlinking is supposed to occur at a very fast time scale, i.e. i>f = e 2 uj and 
v d = e 2 v d , with Uf, u d = 0(1). The macroscopic restoring potential V is defined such 
that: 



Then, 



Similarly, we have 


b(9 1 , 9 2 ) = a | sin(ffi - 9 2 )\ 


a 


e 


sin (0i - 0 2 ) | 
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and consequently, 


dg 1 b(8 1 ,8 2 ) — -do 1 b{9i,9 2 ). 

£ 


Then we have: 


V x , F\ — v^V, 


Xl 


y/eViiVix i, 8 1: £ i, x 2 , 6> 2 , £ 2 )£ 3 5 , (5 i, 0i, 4, ® 2 , #2,4) 


Z£ 


dx^dl x dl 7 


= £ 2 \7x 1 F 1 , 

Flink(x i,0i)= / do 1 V(x 1 ,6 1 J 1 ,X2 : 9 2 ,£2)£ 3 g(Xl,0lJl,X2,9 2 j2) 


dx 2 ^dhdl 2 


J L* c 

. f 1 ~ „ ~ ~ dx 2 —dlidl 2 

F al (x ll 9 1 ) = / -<9 ei 6(6» 1 ,6' 2 )e 3 5f(xi,6» 1 ,£i,x 2 ,6' 2 ,l 7 2 )-^ 

./rx £ £ 

where L e = M 2 x [— |, |] x [— ^|^] 2 . Finally, we define X l and X 2 such that: 


X\(x\, 9i) = y/eX^ixi) + ^£2 S(xi, 0i) = \/£Ad(Xi, 0i), 

£.f 

eddgj + ^eF link ~ 

A 2 (xi, 0i) =-~-= A 2 (xi, 6> 1 ), 

£./ 


with X^A) and X 2 (xiA) defined by 


AAA) = VA°A) + £■j(Xl,0l), 


AAA) = 


ddg 1 f T £F/infc 

~7 ' 


The macroscopic fiber linking/unlinking operator S(g) is similar to the one defined Eq. f!3.5|) . 
Indeed, from Eq. (12.12ft : £(xi, 9 X: x 2 , 9 2 ) = £ _1 / 2 /(xiA, x 2 A) and thus: 

S(g) = F//(xi,0i)/(x 2 ,0 2 )% l!01)5 ; 2ie2) (£)% 1]ei!i2i02) (4) -u d g. 


Altogether, the macroscopic version of Eqs. (14. 1 [) - (14.21) reads (dropping the tildes for the 
sake of clarity): 


A A F al -e ( td 0 F lmk + ddlf ) 




+ £ 2 dtf - V x • (y x Uf) - de(dgUf) - £V X . • F, - dAJ = 0, 


(4.3) 
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and 


~S(g) - £i[d ei (gF a i(x i, 6>i)) + dg 2 (gF al (x 2 ,0 2 ))\ 

—e 2 ^dg 1 (gX 2 (xi, 6 \)) + dg 2 ( gX 2 (x 2 , ^ 2 ))^ 

+£ 3 (d t g-W xl • (gX 1 (x u e 1 ))-V X2 • {gX 1 (x 2 ,e 2 )) 

-dg 1 (gdgU 1 (6 1 )) - dg 2 (gdgU 1 (9 2 )) 

~dX xl • (g^ix^O,)) - dW X2 ■ (g^(x 2 ,9 2 ))j = 0. 

From now on, we note f £ = / and g £ = g. The following proposition holds: 
Proposition 4.1. Assuming f £ and g £ exist, then, formally, they satisfy: 


(4.4) 


-fdg[dg$[f}(x,6)r)-dd 2 e r 


+ £ 


dtf e - • (X x U"f £ ) - d t 


dgU 1 + £G[/ E ](x, 0) 


f £ ~ dA x f £ 


(4.5) 


= 0(£ 2 ), 


and 


g £ ( x 1 , 6 i, l\, x 2 , 0 2 , If) — —f £ (xi, 9±)f £ (x 2 , 02 )de( Xu e u x 2 ,e 2 ) (^)de( X2 , 02 ,xi, 01 )^ 2 ) 

Vd 

+ 0(£ 2 ), 


(4.6) 


with 


<f)[/ e ](xi, 6 >i) = Ci / sin 2 (0 — 9 2 )f £ (x 1 , Of) 


d 0 2 


71 


G[/'](xi,0i)=C 2 £ 

*,J=1 


f £ (% I7 ^2)-Sjj(0i, 02 )-: 

7T 


dxidxj 


Q = £ = 




48 14 * 


and 


-B(0i, ^ 2 ) — sin 2(0i — 02)[ca(0i) (8) ca(0i) + u;(02) <8> uj( 9 2 )\ — [Bij(9\, 02)) ^ y=1 2 - 
Remark 4.1. In the proof of proposition \J~1\ we will show that 

(x, A) = 0(£ 3 ), 

FLkfruOi) = °(£), 

rj(x,,ei) = £ s» 1 t[r](xi.e 1 ) + £ 2 G[/'](x 1 .e 1 ) + o(£ 3 ). 


(4.7) 

(4.8) 

(4.9) 

(4.10) 


(4.11) 

(4.12) 
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The proof of this proposition is given in section 14.31 From these equations, one notes 


that the hypothesis of dominant creation/deletion of links makes the reaction forces F\ 


and Fi ink of order 0(e 3 ). In this case, the process of linking/unlinking is so fast that the 
constraint is satisfied at all times. Moreover, under this assumption, the first contribution 
of the alignment force acting on a fiber is the sum of elementary alignment forces generated 
by its intersecting fibers, weighted by One also notes that the alignment force F a i is 
local in space. 

Under these scaling assumptions, the leading order of the left-hand side of Eq. (j4. 5jl 
takes the form of a collision operator of kinetic theory. It acts on the orientation vector 
6 only and it expresses that the alignment potential (12. 8 p is counter-balanced by the 
diffusion term which tends to spread the particles isotropically on the sphere. The other 
terms act at lower order e. 

As the large scale limit involves an expansion of the solution around a local equilibrium, 
the study of the local equilibria of the collision operator are of key importance. Therefore, 
section [5] will be dedicated to the study of the properties of the left-hand side of (14.5p . 

4.3 Proof of proposition 14.11 

Proof. From Eq. (14.4p . one notes that the source term S(g £ ) is of order 0(e). Thus: 



+ 0(e). 


Inserting this expression into the relations for F[ and Ff ink and Ff d (see Eqs. (13.311 - (13. 4p ). 
one obtains (dropping the tildes for the new variables, and denoting V = V (aq, &i,£i, x 2 , 9 2 , £ 2 ) 
and b e = b £ (9i, 0 2 )): 




Vff £ (x iA) 


A V f £ (x 2 , 6 2 ) 6 i (x , (h) 

(4.14) 


F £ 

r link 





We note that if 0(aq, 61 , £ 1 , x 2 , # 2 , £ 2 ) € L°°(L £ ) with sufficient decay at infinity, then 
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since the measure of L £ intersected with any compact set of M 2 x [—7 t/2, 7t/2] x M 2 is 
of order e. indeed, the domain of integration with respect to £\ or £ 2 has a measure of 
order e. Thus, assuming that the 0(e) remainder in (I4.13|) is an L°° function, which is 
legitimate in view of the diffusive character of (14.4p . we get: 


F[ = ( j V Xl V( Xl , 9 h £(xi, Oi,x 2 , 9 2 ),x 2 , 9 2 , £(x 2 ,9 2 , Xl , 9f) 

rM))<i*^ Y' nxiA) +o&, 

TT J Vd 


F £ = 

link 


d ei V(x 1 , 9 1 J(x u 9 u x 2 , 9 2 ),x 2 , 9 2 , £(x 2 ,9 2 , Xl , 9f) 


K*(x i,0i) 


* e ,_ d0 2 \v f f £ (x 1 ,0 1 ) ^ 

/ (x 2 ,9 2 )dx 2 — 1-h 0(e ), 

7T / V d 


Ki = — f d 6 l b( 0 1 , e 2 )f e (x 2 , e 2 )dx 2 — + 0 (e 2 ), 

v d 7T 

K*(x i,0i) 


(4.15) 


where K £ (x \, 0i) is the set of fibers intersecting fiber in (x, 0), given by: 


K £ (x 1 ,9 1 ) = {(x 2 ,0 2 ) | \£{x 1 ,9 1 ,x 2 ,9 2 )\ < yfeL) 2 , |^(x 2 , 0 2 , ®i, #i)| < yfzL/2}. (4.16) 

From the fact that V is a quadratic function of X\ + l\u>(9f — x 2 — £ 2 co(9 2 ) and the fact 
that setting 0 = £(xi,9i,x 2 ,9 2 ) and £ 2 = I(x 2 ,9 2 ,xi,9i) just cancels this expression, one 
immediately notes that: 


V Xl V(x u 9i,£(xi, 9i,x 2 , 9 2 ),x 2 , 9 2 , £(x 2 , 9 2 , x x , 9f) = 0, 
d 6l V(x 1 ,9 1 ,£(x 1 , 9 1 ,x 2 , 9 2 ),x 2i 9 2 ,I(x 2 ,9 2 , x 1 , 9 1 )) = 0. 


So, finally: 

Ff = 0(e 2 ), FUs,2 = 0(e 2 ). (4.17) 

We are left with: 

F e ai = —f £ (x i,0i) [ de 1 b(9 1 ,9 2 )f £ (x 2 ,9 2 )dx 2 —+ 0(e 2 ). (4.18) 

U d / 7T 

K*(x i,6»i) 


From now on, we write oj\ = uj(9f and u> 2 = oj(9 2 ). By the change of variables x 2 i —> 
(si, s 2 ) defined by 

yfeL yfeL 

x 2 = xi H-—sicni- —s 2 uj 2 , 

with associated Jacobian 

_ Lyfe f cos 9i — cos 0 2 \ 

X2 2 ^siri 9\ —sin 0 2 / 
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and | det(J a;2 )| = ^| sin(#i — 9 2 )\, we have: 


*d(zi,0i) = eC'(a;i,0i) 



sin(6>! - 9 2 )\d 6l b(9 1 ,9 2 ) 


— § |si|,|s2|<l 


/ e (xi + - ^^s 2 uj 2 , 9 2 )ds 1 ds 2 — 

Z Z 71 

+ 0(e 2 ), 


where C^Xi,^) = LzZtLAlllll _ Thanks to (12.71) with jd = 1, one notes that 

<9 0 i&(0i, 0 2 ) = ad 0l | sin(6»i - 0 2 )|, and then, | sin(0i - 6 l 2 )|<9 ei &(6'i, 0 2 ) = f<9 6l sin 2 (6h - d 2 ). 
Then, 


5a 


= -TrC(ziA) 



<%, sin 2 (6 l i - 0 2 ) 


— f |si|,| 82 |<l 


; £ / ^ yfeL Q x , , d 0 2 

J{xi H-sicai-— s 2 uj 2 ,9 2 )dsids 2 - 

2 2 7T 

+ 0 ( 5 2 ). 


(4.19) 


By Taylor expansion, we have: 

f £ {x i + Sicni — ^-^—s 2 uj 2 , 0 2 ) = f £ (xi, 0 2 ) + ^ V x / £ (xi, 0 2 ).(sia;i — s 2 (n 2 ) 

eL 2 

+ (si^i — s 2 o; 2 ) :r V 2 /‘'(xi, 0 2 )(siWi — s 2 cn 2 ) 

+ ^((^ — s 2 (n 2 |) 3 ), 

where V 2 / £ is the spatial-hessian matrix of / e ((V 2 /)ij = g ^J x . ), and for any vector a of 
M 2 and any 2x2 matrix B : a T Ba = 2 j 2 B^ajai. Integrating over s 1; s 2 G [—1,1], 

the odd terms with respect to either Si or s 2 vanish. Therefore: 


7T 

2 

J J d 01 sin 2 (0i - 0 2 )f e (x i + 

“f |S1|,|S2|<1 


\fsL 

— S1W1 - 



S 2 w 2 , 0 2 )ds 1 ds 2 


dd 2 

71 


2 

= 4 f d 9l sin 2 (6 l i - 9 2 )f £ (x i, 0 2 

_ 7T 

2 


f 7 2 7 rW- 

+ — / d dl sin 2 (6»i - 9 2 )X7 2 x f £ (xi,9 2 ) : [cm ® cui + cj 2 ®o; 2 ]— + 0(e 2 ), 
0 7 7T 

_ 7T 

2 


(4.20) 
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where \/A,B e M 2 , A : B = JA je[i 2 ] an d f° r an y vec t° rs £ M 2 , we write 

(a; <S> co')ij = uJiUj'j- Then: 


2 

/ dd 
d 8l sin 2 (6'i - 0 2 )V 2 Xl f e (x 1 , d 2 ) : [cui <g) uq + uq <g) cu 2 ]-^ 

- — 

2 


2 


= E 

hb) =1 


jjE- [ f(x 1 ,e 2 )B ij (91,9,)—, 

C'OCj j OCj J 7T 

_ 7T 

2 


(4.21) 


where: 

Bq(#i,# 2 ) = [^i(^i)wj(^i) + ^ 1 (^ 2 )^(^ 2 )] sin(2(0i — $ 2 )). 

A hrst consequence of what precedes is that = 0(e). Therefore, S 1 ^) = 0(e 2 ) (instead 
of formally 0(e) as seen from Eq. (14.41) ). As a consequence, the remainder in (j4.13[) is 0(e 2 ) 
instead of being 0(e), and the same is true for the remainders in (14.141b Consequently, 
the remainders in (j4.15j) are 0(e 3 ) instead of being 0(e 2 ) as before. It follows that the 
remainders in (j4.17[i - (14.18j) are 0(e 3 ) as well. Then, inserting (14.201) and (14.211) into (14.181) 
(with remainder 0(e 2 )), we get (j4.6j) - (|4.5j) . which ends the proof. □ 

From now on, we focus on Eq. (14. 5 j) in which we neglect the 0(e 2 ) terms, namely 


-adel W(i,e)/ e ) -ddlr 


+ £ 


dtf e - v, • (V x U°n - d ( 


deU' + ZGirj&e) 


r - dA x r 


= 0 , 

(4.22) 


where <f> and G are given by (14.71) - (14.81) respectively, and we investigate the limit e —> 0. 
This is the object of the next section. 


5 Large scale limit 

In this section, the limit e — > 0 of the solution f e to (14.22ft is explored. For this purpose, 
Eq. (I4.22j) is rewritten 

dtf E - v,. • (W x U°n - d e ((d e U 1 + £G{f])f £ ) - dA x f = (5.1) 

where the collision operator Q(f £ ) is defined by 

Q(f) = dd 2 e f + Zd e (d e m)f ), ( 5 - 2 ) 

r f jfl 

®lf} = Ci sin 2 (0 - d 2 )f— , (5.3) 

J -E 7T 
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and where we recall that C\ and G[f] are dehned by (I4.8[) and (14.91) respectively. The 
operator Q is a non linear operator on / which acts on 6 only and leaves x and t as 
parameters. For each function <E>(0), we define M§(6) by: 

Af*(0) = i e -^ (0)/d , (5.4) 

Zj 

7T 

where Z is a normalization factor such that Z = fJ E e~^ e ^ d ^-. Thus, M$(0) is a 
probability distribution of 6. Such functions are called generalized Von Mises distributions 
(the Von Mises distribution being the case of $(0) = — cos 9). The next section is devoted 
to the analysis of the properties of Q(f) and follows closely Ref. [IT3] . 


5.1 Properties of Q 

5.1.1 Equilibria 

In this section, the equilibria of the operator Q are studied, and the following proposition 
is proven: 

Proposition 5.1. Here, we restrict ourselves to functions of 9 only. 

(i) The operator Q can be written: 

QU) = dd e (M m d e (jf—j)- ( 5 - 5 ) 

(ii) The equilibrium solutions of Q, i.e. the functions f such that Q(f) = 0 are of the 

form f{6) = where is defined by Eq. (j5.4[) and p is a positive constant. 

This proposition shows that the equilibria of operator Q are generalized Von Mises 
distributions of 9, weighted by the particle density. 

Proof. To prove (i), one can note that: 

dd e = dd e (W - fd e (\og(M m ))^ = d ( > (dd e f + £<%$[/]/ 

= W)- 


To prove (ii), note that / = is solution of (15.51) . Conversely, suppose that / is such 

that 

^ a 4 ) ) =o - 

We define the sets Hf and Vf by: 


Hf = {c f> measurable on [— 


7T 7T 
2 ’ 2 



\ r d9 
M $\/]— < +oo}, 


do 


'M, 


* 1/1 


d9 

M$[f]— < +ooj. 
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and 


V f = {(j) e H 



















The norms 


| H f , || • || v f 011 Hf and Vf are then defined such that: 
Mv, = + \<l>\v f - 


where 


and 


11010 / - 


rir/2 

0 

J —7T/2 

M $[/ ] 


/•tt/2 
/—7T/2 


da 


101 V> = 

For / G Vf using Green’s formula, we get: 

/ dd 


'M, 


r 7r / 2 / f 

ddgi M m d e {- 


I —7r/2 


'M, 


m 


M m * 


m 


= -d 


dd 

M m—> 


dO 


f 71-/2 


M, 


' —7r/2 


m 


mJ—) 


'M. 


*[/] 


'do 


7r 


= 0 , 


and thus, dg( ) = 0. Then, / = with p > 0, which ends the proof. 


□ 


Now, the following lemma is proven: 

Lemma 5.2. For any function f(0), the potential function <E>[/](0) of Eq. (I5.3jl can be 
written: 


Wm = C- < F Vf cos2(6-e f ), 
r 71- / 2 r dd 


(5.6) 


where C\ is given by (14.91) . C = C - J p-, Pf = fl^f^F an d 0/,0/) £ x [—f, f) are 
uniquely defined by: 


Vf 


cos 2 Of \ 
sin 2 Of 


' 2 (cos 20'\ dO’ 

, sin 28' ) J[l ’ vr : 


or equivalently by: 


‘2 A A' 

COS2 = 

_ 7T 7T 

2 


a nr 

sm2(0' -6 f )f(0')— = 0. (5.7) 

7T 


Remark that the second condition is equivalent to saying that 

1 _ w /sin207(0 , )d0' 

2 taT 


0j = - tan " 1 ( J r 0111 

J O V f oa/ r QI\UQI / ’ 


/cos 20'f(0 , )d0' 

and this defines Of uniquely modulo n. 

Proof. As sin 2 (6 l — 0') = |(1 — cos 20 cos20' — sin 20 sin 20'), <£>[/] can be decomposed into: 

d& 


*m) = Ci sin 2 (0 - 0')f(0') — 

I 7r 7r 


C/ 

2 


00' 


dO ' 


/(0')-cos20 / cos20'/(0 / )-sin20 / sin20 / /(0 / ) 


dO’ 


7r 


7T 


7T 


= y(p-0 /C °s20-^/)), 


The result follows. 


□ 
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Let us now suppose that ^ depends on rjf: 

Hypothesis 5.1. The parameter ^ is supposed to be inversely proportional to the local 
fiber density: — = —, with 7 a constant. 

J * 1 /, nt 7 1 


Note that, thanks to Hypothesis 15.11 we have 

jCirjf _ jL 2 v f 1 _ £oL 2 7 _ 

2d ~ a 2u d Vf 2d~ 4d ~ ^ 

where r is a constant depending only on the data of the problem. 


(5.8) 


Proposition 5.3. Here, we restrict ourselves to functions of 9 only. Under Hypothe- 
sis 1 5.1\ the equilibrium solutions of Q, i.e. the functions f eq such that Q(f eq ) = 0 are of 
the form: 

f eq (9) = pM do (9), (5.9) 

for arbitrary p G [0, 00 ) and 0 o G [—|, |) and where: 


( 




3 r cos 2 ( 0 — 9q) 


Me 0 - 


Z 

Z = Z(r) = 




c r cos 2 ( 0 - 0 o ) d® 
TT 


(5.10) 


with r given by (I5.8[) . We have rjf = pc(r) with 


c(r) 


f 2 n cos 26e 

J 2 


r cos 26 d6_ 

7T 



gr cos 26 d§_ 

7T 


(5.11) 


Proposition 15.31 gives a precise description of the equilibria of Q, in terms of classical von 
Mises-Fisher distributions. 


Proof of proposition Id.51 From Proposition 15.11 the equilibria of the collision operator 
Q(f) are of the form 


f = P 


= 5 d 


f 


f 2 n e 

J 2 


Thanks to Eq. (14.91) . Lemma l522l Eqs. (15.8ft and (15. lip , we get: 


m=p 


cos 2(0 6 f ) 

f\ e -^T+^rVfCOs2(8'-e f )dW 
2 n 


= p(x) 


g r cos 2(6—6 f) 


f 2 £r cos 2(6'—6 f) d6' 

■J — 77 7T 


(5.12) 


where {rjf, Of) G M + x [—§,§) satisfy Eq. (15.71) . Therefore, / is of the form (15.9j) with 
r = e,C ^ s . By Hypothesis 15.11 and (14. 9ft . r = ■ Conversely, let / be given by (15.91) . 

Then, by (I5.6p .and (I5.8j) . 4>[f] = C — r| cos 2(0 — Of) with Of uniquely determined by 

/^ 2 sin 2(0 - 6f)f(6)f = 0. But ff( 2 j2 sin 2(0 - 0 o )f(0)f = 0 by symmetry, showing 
that Of = 0 o mod(7r). Therefore, M^p = Mg 0 and / = pM^p showing that / is an 
equilibrium, which ends the proof. □ 
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Thanks to Eq. (15. lip . Hypothesis 15.11 amounts to supposing that the ratio — is in- 

^d 

versely proportional to the fiber density. 

Since there is no obvious conservation relation other than the conservation of the local 
fiber density, the only collision invariants in this model are the constants. The integration 
of equation (14. 5 p against these invariants does not allow us to find the evolution equation 
for the mean orientation. In order to obtain an equation on 9q, inspired from Ref. 53, the 
concept of Generalized Collision Invariants (GCI), i.e. of collision invariants when acting 
on a restricted subset of functions /, is introduced. 


5.1.2 Collision invariant 


A collision invariant is a function T such that for all function / of 9, f Q(f)'i>d9 = 0. 
However, due to the lack of momentum conservation, the only collision invariants are 
the constants. This is not enough to determine both p and 9 0 . To this aim, following 
Refs. [T9] and [T7], we introduce the notion of GCI. For any 9 0 G [—f1), we define Lg 0 as 
the following linear operator: 


Le 0 f — dd g [ M eo d g (Y—) 

Mg 0 


Note that Q(f ) = Lg f f where 9f satisfies Eq. (I5.7p . 

Definition 5.4. For a given 9 0 G [— f,f) a GCI associated to 9 0 is a function T such 
that: 

/ Lg 0 f T— = 0 Vf such that 9 f = 9 0 mod(7r). (5.13) 

/ vr 7T 

J ~2 

The set of the GCI associated to a given 9 0 e [—|, |) is a linear space denoted by Qg 0 . 
Lemma 5.5. T G Qg 0 if and only if 3/3 G R such that: 

L* do y = Psm2(9-9 0 ), (5.14) 

where L* e is the L 2 formal adjoint of Lg 0 , i. e. 


d 


L* 0 T = 

e ° Mf 


e 0 


do Mg 0 dg T 


Proof. By (15.7p . the condition 9f = 9 0 mod(7r) is equivalent to the linear constraint: 

/•f d9 

/ / sin 2(0 — 9 0 )— = 0. 

J-l 71 

By a classical duality argument 53. we deduce that T G Qg 0 if and only if: 


3(3 G M such that 


d9 


d.9 


Lgjy- = f3 / / sin 2(9 — 9 0 )— V/. 


IT 


IT 


Note that now, there are no more constraints on /. Therefore, we can eliminate / and 
get tfEUD. □ 
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Proposition 5.6. Any GCI ^fg 0 associated to 9 0 can be written: 


^e 0 (9) — c + ~ Oo), 


(5.15) 


with arbitrary C, /3 £ 
whose expression is: 


R and with g an odd tt periodic function belonging to Hq(0, |) 7 


m 


1 

2 r 




—r cos 2 9' d,0' 



26' d6' 


(5.16) 


Proof. Following Refs. [19], [IT], using Lax-Milgram’s theorem and Poincare’s inequality, 
it is easy to show that the problem L^T) = |/3 sin 2(9 — 0 O ) has a unique solution in the 
space H 1 ( — |, of functions H L (—^, |) with zero mean. Then, the change of variables 
9' = 6 — 6*o is performed, and functions of the form T(0) = (3g(9) with g odd are searched. 
Then, T £ FT 1 ([— |, |]) if and only if g belongs to 77 q( 0, |). Straightforward computations 
show that T is a solution of (15. 14(1 if and only if g is a solution of 


(M 0 g')' = -sin20M o . (5.17) 

As M 0 (9) = erc ° s2a and as we search for g £ 7/q( 0, |), an analytic expression for g can 
be found. Indeed, since — sin 2 9M 0 = f^M 0 , integrating 05.1 7p with respect to 9 once, we 
get: 

g'(9 ) = — + CZe~ rcos2e , 

2 r 

for an appropriate constant C. Then, since g £ T7o(0, §), 


9(«) 


2f + CZ 


—r cos 20' 


d0'. 


Finally, as g £ 77g (0,7r), g(0) = g(7r) = 0 and C can be determined: 


C 


71 

ArZ f ( f e~ rcos26 ' d9' 


1 

2 rZ f^ 2 e- rcos20 'd9' 


1 

2rZ 2 ’ 


Indeed, we have: 


t/2 


7T/2 


—r cos 20' j/V 


d0' = 


t/2 


'-7T/2 


cos 20' j/V 


d0', 


by the change of variable 0—* | — 0 for 0 > 0 and 9 —* — 0 for 0 < 0. This yields the 

result. For further usage, we note that 


m 



(5.18) 

□ 


25 












5.2 Limit £ 


0 


In this section, the formal limit e —> 0 of Eq. (14.51) is studied. We aim to prove the 
following theorem: 

Theorem 5.7. Under the scaling \472 \ and \5.1[ the solution f £ of eq. (15.1ft formally con¬ 
verges to f(x,9,t) given by 

f(x,0,t) = p(x,t)M eo [ Xft )(9), (5.19) 

where M 8q is given by (15. 10[) and p(x,t) and 9 0 (x,t ) satisfy the following system: 

dtp - V, ■ (' V x U°p ) - dA x p = 0, (5.20) 

and 

pdtOo ~ pV X U° ■ VA — 2a 2 V x p ■ VA — a 2 pA x 9o 

+a 3 (p\7 2 x 9 0 + VA 0 W x p + V x p 0 VA) : [u 0 0 cu 0 - Uq 0 A] (5.21) 

+ (2pa 3 VA 0 VA - aqV^p) : A 0 w o 1 + w ( (® A + a b p(d e U 1 ) = 0 , 

7t/2 

where ( h) = f h(9)Mg 0 (9 )^ for any function h of 9 e [—§,f), and where the coeffi- 

—7r/2 

cients a 2 , a 3 , aq, aq are given by: 

d t , £oL 4 7c(r) 

“ 2 = A + 24d >’ 

_ £aT 4 7 ( 1 , 6 dc(r)^ 

a 3 — -ttt:—(777J — 1 + — r9_ )j 


oq — 


24a x HZ 2 
^aL 4 7 


faL 2 ^' 


with 07 given 6 g: 


192Z 2 Q! 1 ’ 
1 

«5 — -j 


a± = 1 


(5.22) 


1 

Z 2 ' 


(5.23) 


Proof. Suppose that all the functions are as smooth as needed and that all convergences 
are as strong as needed. In the limit e —> 0 , let f £ —$■ f. As Q(f £ ) = 0(e), then Q(f) = 0 . 
By proposition 15.61 we deduce that / is given by (15.19j) with p > 0 and 9 0 e [— f) to 
be determined. In order to find the equations for p and 9 0 , we use the set of GCI given 
by Prop. 15.31 


Equation for p The use of the constant GCI amounts to integrating Eq. (15.11) over 
[-f, §). This gives: 


dtf e - V x • (VA7 e ) - d e 


dgU 1 + fG[f £ }(9) 


f ~ dAJ 


r , d9 
£( - = 0 , 
7T 


which leads to the continuity equation for p £ \ 

d t p £ - V x • (' S7 x U°p £ ) - dA xP £ = 0. 
In the limit e —> 0, p £ —> p which leads to Eq. (15.20(1 . 
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Equation for 0 O We multiply Eq. (15.11) by the GCI '•l ) o f£ associated with the direction 
dfe of f £ , namely = p(0 ~ 0/®) where g is the function defined in Prop. 15.61 We 
integrate with respect to 6 and first note that: 

7T 7T 

2 Q(f e )*o fe d6 = 1 2 Le f J £ ^e S edO = 0 , 

/ 7T / 7T 

J ~2 J ~2 

by (I5.13p . Since f £ —> pM 6o , we have 0/e —$■ 0 O and 4d /e —$■ d/ 0Q . Therefore, in the limit 
£ —y 0, we get: 


d t (pMe 0 )-V x iV x U 0 pM eo )-d e 


dgU'+tGlpMeM 


pM 0o ) - dA x (pMg 0 ) ) T 6o d0 = 0. 

(5.24) 


For simplicity, we denote Mg 0 = M. We have: 

A x (pM) = MA x p + pA x M + 2V x p • V X M, 

V x • C V x U°pM ) = MV, • (V x U°p) + pV x f/° • V,M. 

Using the continuity equation fl5.2()[) . we have: 

d t (pM) = pd t M + Md t p = pd t M + (V, • (V X UV) + dA x p)M. 

So: 

dt( P M) - V x • (' V x U°pM ) - dA x ( P M) = pd t M - p\7 x U° • V X M - dpA x M - 2dV x p • V X M. 
Therefore, Eq. (15.241) reads: 


p [ 2 — - Ad - X 2 - X 3 - Ad = 0, 


where: 


Ad = 


, V x p, 


p(V x U u + 2d——) • V X M T 

P / 7T 


d6 


do 


x 2 = I "deldoU'pMjV—, 

X3 = ^f] de ^ 

/■f dd 

Ad = dp / A X MT —. 

/ 7T 7T 

J ~2 


We now turn to the development of each term of Eq. (I5.25p . We have: 

V X M = 2 r sin 2(0 — 0 o )MV x 0o. 

Then, 

(V X U° + 2d—• V X M = 2r sin 2(0 - 0 O )M( V X U° + 2d— ) ■ V X 0 O , 
P \ P J 


(5.25) 

(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 
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and thus, X\ can be written: 


From integration by parts, the following relations can be written: 

< S in2(0-(W*) = X (1 _J_ ) = _L Q , (5.31) 

Therefore, we have: 

x, = ^(X x U° • V x 9 0 + 2 d WxP ' Wx9 ° ). (5.32) 

2r p 

Since X 2 is the integral of a 7r-periodic function over a period, we can write 

/■<9 0 +7r/2 r n 

X 2 = / do&U'pM)^ — . 

J 9q—-k/2 71 

Now, by construction, (see prop 15. 6p . T(6 I 0 — |) = \P(0 O ) = T(6 l o + |) = 0. So, integrating 
by parts, we have 

f0o +f c m 

X 2 = - pM deU 1 <%T —. 

J Qo—tt/2 ^ 

Now, by construction again (see (I5.18P ). we have 


d e q = — (1 - 1 -). 

2 r V MZ 2 ’ 


Using again the 7 r-periodicity of U 1 , we obtain: 

p 1 


X 2 = -^-(dgU\ 1 

2 r 


MZ Z = -> U ' 7 


(5.33) 


(5.34) 


Now, let us turn to X 3 . The details of this computation are postponed to appendix] 
We find: 


X 3 


dL 2 

~\2 


- c(r)(pA x 9 0 + 2V X 6» 0 • V x p) 

+ (2^71 V x 6> 0 0 V x 9 0 - -^V x p) : [wq (8 ) 0 w 0 ] 

+ 7i(pV x 6> 0 + V x 9 0 0 V x p + V x p 0 V x 9 0 ) : [ca 0 0 u 0 - 0 a;;}] , 


(5.35) 


1 3c(r) 1 6 dc(r) 

4Z 2 ~ + 2r ~ 4Z 2 ~ + aL 2 i 7 ' 


where, using (15.81) . 


















r 4 -i 

We note that ^P -71 = j-a 3 . Finally, let us explicit the last term X 4 . A direct computa¬ 
tion gives: 


A X M = M 


4r[rsin 2 2(0 — 9 0 ) — cos2(9 — 0 O )] I V^o ) 2 + 2r sin 2(0 — 9 0 )A X 9 0 


Then, we deduce that 


X 4 = dp 


2 rA a , 6 , 0 (sin 2(9 — 6 ) 0 )'F) 


+ |Va;0o| 2 4r[— (cos 2(0 — 0 Q )\k) + r(sin 2 2(9 — 6 *o)T)] 


By symmetry, we have: 

(sin 2 2 ( 0 - 0 o )T) = i(cos 2 ( 0 - 0 o )T). 
r 

Therefore, with (15.311) . we get: 

X 4 = ^(1 - ~)A X 9 0 = ^AJo. (5.36) 

Now, d t M = 2r sin 2(9 — 9 0 )Md t 9 0 , and 

7r 

p I d t M T = 2rp(sin2(6> - 9 0 )'£>)d t 9 0 = ^-(1 - -^)d t 9 0 = ^d t 9 0 (5.37) 

Collecting (I5.32f) to (I5.36P and inserting them into (15.25j) leads to (I5.2ip . □ 


6 Case of a homogeneous fiber distribution: station¬ 
ary solutions 

In this section, we study the stationary solutions of (I5.20p - (l5.2ip in the case of a spatially 
homogeneous fiber distribution and consequently no external spatial potential U° — 0. 
We make the following assumption: 

Hypothesis 6.1. The fiber spatial distribution is supposed to be homogeneous, i.e. there 
exists a constant p 0 > 0 such that p(x, t ) = po for all (x, t) G R 2 x [0, 00 ). We also suppose 
that there are no external spatial forces, i.e. U° = 0. 

We first note that in the absence of external forces, a uniform and constant density 
Po is a solution of Eq. (I5.20p . Now, we are interested in the stationary solutions for the 
fiber orientation equation (15.211) . Noting that the terms involving the spatial derivatives 
of p, we End that such stationary solutions satisfy the following equation: 

oi2A x 9Q—afiuQ ( 8 ) cuo — uyj - ® : V 2 $o 

( 6 . 1 ) 

- 2 a 3 [oj 0 <g> uj 0 +uj 0 ® £u 0 ] : V X 0 O <g> \7 X 9 0 = a 5 (d g U ). 
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In this equation, the coefficients r, aq, a 2 and 03 are constants thanks to (15.81) . Moreover, 
using (I5.22[) . they can be written as functions of d, L 2 and r as follows: 


011 (r) 
< a 2 (d,r,L 2 ) 
a 3 (d, r, L 2 ) 


d\ 1 + 


Z(r) 2 ’ 
L 2 rc(r ) 


dL 2 r 
6 aq (r) 


6 aq(r) 

A(r). 


( 6 . 2 ) 

(6.3) 

(6.4) 


with 


A(r) 




(6.5) 


We now show that (16.ip is an elliptic equation. We first introduce some definitions. 

Given a function f(x, E ) smooth in its arguments x G hi, £ 6 1 x I ] x S 2 (R), where 
5 2 (M) is the space of 2 x 2 symmetric matrices with real coefficients, we dehne the non 
linear differential operator F : C'°°(M 2 ) —y C' 00 ^ 2 ) such that for any x G K 2 and any 
u e C°°(K 2 ), we have 

F(u(x)) = f (x, D 2 u(x)), 


where D 2 u = {D a u,a e N 2 , |a| < 2} and where, for a multi-index a = (ai,a 2 ) G ff 2 , 
|a| = aq + a 2 and D a u = jSrj^k ■ The operator F is said to be elliptic at u\ G C°°(K 2 ) 

(see Ref. [34]) if its linearization DF{u \) is an elliptic, linear differential operator. We 
state the following proposition: 


Proposition 6.1. Eq. \6.1\) can be put in the form 


f(x, D 2 9 0 (x)) = 0, x G K 2 , ( 6 . 6 ) 

where f(x, D 2 9 0 ) is the following operator, quasi linear in 9 0 : 

2 

f(x, D 2 9q) = d Xi (a ij (9o)d Xj 9o) - a 5 h{0 0 ). (6.7) 

ij =1 


Here, h(9 0 ) = ( dgU 1 ) and A{6q) = (ajj(# 0 ))ij=1,2 is a 2x2 matrix such that: 




—ck 3 sm26 , 0 


Moreover, if the following condition is satisfied for all r G 

A(r) + c(r) > 0, 


( 6 . 8 ) 


(6.9) 


where A(r) is given by Id. 51) . then F{9) = f(x, Zl 2 ^) is elliptic at 9i for all 9\ G G 2 (M 2 ). 
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Proof. For any # G [— |,~), letting u>(9) = (cos#, sin#) and a;- L (#) = (—sin#, cos#), we 
have: 

■^-[u(9) <g) ca(#) - o;- L (#) <g> cj ± (#)] = 2[u(9) <g) cj- l (#) + w ± (#) <2> u(9)\. 
du 

Let F : C' 00 (IR. 2 ) —>■ C^R 2 ) be the non linear differential operator defined by: 

F(# 0 ) = f(x, J D 2 # 0 ), 


for f defined by (I6.7|h Let DF(6i) denote its linearization at 6 \. Then, DF(9i) is a linear 
map from C 2 (R 2 ) to (^(R 2 ) and reads, for v G C 2 (R 2 ): 


DF(6 i)v 


dF(9 1 + sv ) 


ds 


s =0 



+ Lv, 


( 6 . 10 ) 


where L is a linear differential operator of order 1 the coefficients of which depend on 
D9\\ 


L v = (^'ij^i){d Xi 9 1 d Xj v + d Xi vd Xj 6i) j + ^ 9 Xi U°d Xi v 

i,j=l ' ' i=1 

+ (a''(#i)^#i^#i + a! i: j(0i)d XiXj 6i — a§h\0i)\v , 

m=i ^ 2 

where aL(#i) and a"(#i) are the first and second order derivatives of the coefficients of 
matrix A which read: 


( a p'(#l))jj = l,2 

< 


/ sin 2 #! — cos 2 #!\ 

y— cos 2 #i — sin 2 #i J 


f cos 2 #i sin 2 #! \ 
y sin 2 #i — cos 2 #iy 


5 


Therefore, the linearization of F at #i is elliptic provided that the matrix Ll(#i) = 
KW)),,,.!. 2 is positive-dehnite. 

Note that the determinant of the matrix A(# x ) = is given by det(A(#x)) = 

«2 — and does not depend on #i. Moreover, det(A(#i)) > 0 provided that |^|| > 1. 
The eigenvalues of the matrix A(#i) solve 


det(Ll(#i) — XI) = A 2 — 2 A «2 + ct 2 — cc 2 = 0 


and the determinant A = 4a 2 is strictly positive as long as a 3 A 0- I 11 this case, the 
matrix A(#i) has two distinct real eigenvalues given by: 

A ± = («2 ± 03 )- 


Therefore, the matrix 7L(#i) is positive definite if and only if 0:2 > |«3 
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We now analyse the sign of each coefficient 01 , 0 : 2 , 0 : 3 . First of all (see Eq. (15. lOjl ). the 
/nth derivative Z^ of Z with respect to r reads: 

7r/2 

Z {p \r)= /(cos^e—f, 

—7T /2 

and we have Z^ 2 k 2 > ( r ) > 0 for all k G N + and all r G M + as the functions 0 —>■ (cos 29) 2k e rcos2e 
are positive for any r G M + . We deduce that Z^ 2 k+ 1 >{r) are increasing functions of r for 
any k G N + . Note that from the symmetry of the function cos 26, we have for any k G N + : 

7r/2 

J (cos 2 £) 2fc+1 ^ = Z ( 2 fc+ 1 ) ( 0 ) = 0 . 

— 7T/2 

Therefore, we also have that ,gl 2 fc+ 1 )(r) > Z^ 2 fc+ 1 )(0) > 0 for any k G N. We thus obtain 
that for any p G N and any r G M + : 

Z (p) (r) > Z {p \ 0 ) > 0 , 


and we note that Z(r) —* cx) as r —* oo. Moreover, as Z( 0) = 1 we deduce Z(r) > 1 for 
any r G [0, +oo). We also note that: 


c(r) 


£<%) zWp) 

Z(r) ~ Z(r) - ’ 


and we have: 


ai(r) > 0, a 2 (d, r, L 2 ) > 0 V(r, L, d) G M + x M x 
Now, by integration by parts, we can write: 


c(r) 
r rZ(r ) 


r/2 


cos 26*e 


r cos 20 


do 


t/2 


2 r%n^r cos 20 


— 7r/2 


7 T Z(r) 


sin 2 20 e 


dO 


IT 


— 7r/2 


t/2 


= i - 


2 cos 20 


Z(r) 


cos 26e 


d6 


= 1 - 


7T 


— 7t/2 


Z <2) (i 


We now show that 

or, equivalently, that 
Indeed, 


Z^ 2 \r) ^ 1 
Z(r) ~ 2’ 

Z(r) < 2Z (2) (r). 


d6 


dd 


d6 


Z(r) — e rcos2 ”—= cos 2 26e rcos20 — + sin 2 26e rcos20 —, 

J TT I 7r ./ n 


( 6 . 11 ) 


( 6 . 12 ) 
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and, by integration by parts, we have: 

7T 7T 7T 

2 2 2 

[ sin 2 20e rcos20 — = / cos 2 20e rcos2e —- r [ sin 2 20 cos 20e r cos20 —. 

./ 7T ./ 7T ./ 7T 


To show that J sin 2 20cos20e rcos20 — is positive, we can note that it is an increasing 

_ 7T 

2 

7T 

2 

function of r and that for r = 0 we have f sin 2 20 cos 20— = 0. Indeed, the derivative of 

J 7T 

_ 7T 

2 

this term with respect to r reads: 

7T 7T 

2 2 

— ( [ sin 2 20 cos20e rcos20 —) = / sin 2 20 cos 2 20e rcos20 —, 
dr v ,/ 7 t J 7r 

_ 7T 7T 

2 2 

which is positive for any r > 0. Therefore: 


,00 


d9 


sin 2 20e r cos 20 — < / cos 2 20e r cos 20 — 

7r ./ 7T 


for any r > 0, and inserting this expression into Eq. (16.12[) . we obtain: 

7T 

,d9 


Z(r)< 2 I cos 2 20e rcos20 — = 2Z (2) (r). 

7T 


All together, we have: 


c(r) 1 
r - 2’ 


for any r £ M + . This relation together with the fact that 1 — j^hyj — I ^ ea ds to: 


, r2 . dL 2 r ,3c(r) . 

a 3 (d,r,L 2 ) = - (1 


6«i(r) 2r 


4Z 2 


))<o. 


Now, we can write: 
or equivalently, using Eqs. 


I ^ I 

CC3 

and 


>1 -v=> o?2 > — CC3, 


i a 2 , ^ 6«i(r) ... 

W >: w > - Wr)+c(r)) - 


Therefore, if (16.9[) holds, then cc 2 > 1 0 : 3 1 and the matrix A(9i) is positive definite for all 
r £ [0,1], L £ R + , d £ M, independently of 0i £ C 2 (M 2 ). We conclude that F is elliptic 
at 9\ for all 0i £ C 2 (M 2 ), provided (16. 9 b holds. □ 

Remark 6.1. As shown by Fig. [H ; A(r) + c(r) is positive for any r £ R + . The rigorous 
proof of this fact will be the subject of future work. 
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Figure 2: Functions A(r) (black), c(r) (green) and A(r) + c(r) (red) as functions of 
r G [0,100], 

7 Conclusion 

In this paper, we have formally derived a macroscopic model for temporarily linked fibers 
interacting through alignment at the links. We have shown that the corresponding kinetic 
model involves two distribution functions: the fiber distribution function and the cross¬ 
link distribution function. The latter can be seen as a joint two-particle fiber distribution 
function. This model provides a unique explicit example of a kinetic model closed at 
the level of the two particle distribution function. We then considered the regime of a 
fast fiber linking/unlinking process, where the link distribution function can be expressed 
simply in terms of the fiber distribution function. We studied the diffusive limit of the 
resulting equation and obtained a system of two coupled nonlinear diffusion equations for 
the fiber density and mean orientation. In the homogeneous fiber density case, we showed 
that the resulting quasilinear problem is elliptic. Future works will deeper investigate the 
mathematical properties of the models, such as rigorously proving the mean-held kinetic 
limit of the particle system or proving existence and uniqueness of smooth solutions for 
the macroscopic diffusion system. Numerical simulations will be performed to validate the 
macroscopic model by comparison with the individual based model. Further perspectives 
are the removal of the fast fiber linking/unlinking hypothesis, in order to understand how 
a finite lifetime of the cross-links affects the macroscopic dynamics. 
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A Proof of Theorem 13.1 


A.l Evolution equation for the fibers 

For all observable functions <t>(x, 9), we define: 


N 




i —1 


Similarly, for all two-particle observable functions ^(xi, 9 1 , £±, x 2 , 9 2 , £ 2 ), we define: 

((g K ,^)) = I ^(x 1 ,6 1 J 1 ,X2,92,£ 2 )g K {x 1} 9 1 ,£ 1 ,X2,92,£2)dx 1 dx2——d£ 1 d£2 

J vr 7 r 


1 

2W 


K 


k =1 




+ 0j(k)> g<k)’ Xi(k), 0i( fc)> Ofc)) ) > 


where integrals over x are carried over M 2 , in 9 over (—|,f) and in £ over ( — §, §)• 
We recall the notations <7j fc)j . (fc) = (X i(fc) , % fe) , ^ (fc) , X m , 9 m ,£ k j{k) ) (resp. C^ k)m = 


(Xj(k)i dj(k)i£j(k)iXi(k), 9i(k),£ i ( k ' ) ))- Then: 


d 

dt 


(/", 4) = ^ £ (vMXitt ),0i(t)) ■ + deHMt), 0M)^r 

1=1 ^ 


Using (12.1 Op and (12. lip , we obtain: 

d 


dt 

1 

'n 


<r» 


1 N 

-Y 

AT 


i= 1 L 


(/iV x $ ■ V,U + Ac«f/)(AA, 0*) 


+d(//V*$ • V ;c log(/ 7V ) + XdgQde log(/ JV ))(W, Oi) 

1 K 

+gVMXi,9 l ) ■-YVx.VdmW + V, 2 U^ (fe) (f))(^ fe)jj(fe) ) 


k =1 






fc=l 


if 


Ac^t&pQ, 6*j)- di(k)(i) + dg 2 b ^j(fc) (*))($»(*;), ^j(fc)) 


fe=i 
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We get, using the definition of a distributional derivative: 


4 </ k , j >) = ( t iV x -(f rt v x {u + d\o g f N )) + xa e (f N a e (u + d\ogf N )),t) 


~Tn £ ("<)' £(v»ri i(t ,(i) + v n v5 m mcm, m ) 


N 


2N 


i =1 


ic 


fc=i 


A 


N 


i =1 


A 

2 N 


K 




fc=l 






£ a.4>K,« i )£(9.,<> 


i =1 v fc=l 

Now, exchanging the sums in i and k in the previous equation, one obtains: 
d 


dt 


(f N ,$) = (ii\7 x -(f N V x (U + dlogf N )) + \d e (f N de(U + dlogf N )),<S>) 


/i 


K 


2N 


k =1 


E v xi v(c? (kim ) ■ \/.Mx m ,9 m ) + w X2 v(c^ m ) ■ vMx m: e m )) 


K 


A 


k =1 
K 


— E ( ^^(C* (fc)J(fc) )9^(X i(fc) , 0 i(fc) ) + de 2 V(C* {kim )d e <S>(X m , 0 m )) 


"y 1 ( ®0ib(@i(k)i 9j(k))deQ(Xi(k), $*(&)) + de 2 b{9i{k)i @j(k))deQ(Xj(k), @j(k)) 


2N 

k=i 

From the symmetry of V (see Eq. (12.41) 1. the following expressions hold: 

v„n<W)) = d».v(c^ m ) = d ei v(c" mm ), 

and from the symmetry of b, we have: 

de 2 b( 0 i(k), 9j(k)) <9$!), 

leading to: 

= R-(fV,((/ + d log/ v j) + \de(f N d„(U + d log/")),*) 


w£ ( v » y ( c hww)' W4>(V,+ V„V{C% mk) ) ■ VMX m ,» m )) 


K 


2N 


k =1 




fc=i 


-^£ ( 9 <-V(C'5 t)JW ) 9»4>(V W ,9. W ) +ftV ( Cyi, i(t ,)a„1> ( A', w ,9, w )) 


A 




fc=i 
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or again: 


J t (f N ^) = (^ x -(f N \7 x (U + d\ogf N ))+\de{f N d e (U + d\ogf N )),<S>) 


K „ 1 


K 


“jvfe S + 0 1,-^2, 02,^2), 


k =1 


(/^V xl V(xi, #i, X 2) 6*2, 4 ) ' V^XiA) 

TAd^ V(xi, 6>i, £i, x 2 ,6*2, £ 2 )d 8l ^(x i, 0i) 

+ A9 ei 6(6»i,6'2)(9 ei $(a:i,6'i)))). 


Therefore, we obtain: 


dt 


(f N ,$) = (^ x -(f N V x (U + dlogf N )) + \d e (f N d e (U + d\ogf N )),<S>) 
-jj{{g K ,^x 1 V(x ll 0 1 ,£ 1 ,X2,9 2 ,£2) ■ V x $(xi, 6i))) 

— ~jy((g K i A r^ 1 V( 2 :i, &i, £i, x 2 ,9 2 , £ 2 ) + d gi b(6i, $ 2 )^ de*&(%h 0i)))- 


Finally, we get: 

j t (f N ,$) =(^ x -{f N V x (U + d\ogf N ))+\d e (f N d e (U + d\ogf N )),<f>) 

+ ^((V xl -(g K \7 xl V),$(x 1 ,9 1 ))) 

+ ( g K de 1 V + d 9l b), $(xi, 0i))). 

=(^ x -(f N V x (U + dlogf N )) + \d e (f N d e (U + dlogf N )) 

+ Jf[[^xi ■ (g K V xl V) + A d 0l (g K {d 9l V + V))]]^, 0i), $(xi, 0i)), 

(A.l) 

where, for a distribution T acting on functions of (xi,9i,£i,X2 } 9 2 ,£ 2 ), we denote by 
[[T]](xi, 9i) the distribution which to any function <&(a;i,0i) associates 

{[[r]](i 1 ,e 1 ),<i>y 1 ,() 1 )) = {{T,ti)), 

and where 1 is the constant function of the variables (aq, 9i, £ 1 , x 2 , 9 2 , £ 2 ) equal to 1. In 
the formal limit N —> 00 , jj —> £ and given the assumptions on the regularizing sequences 
£ N , V N , we get that f N ->• /, f N ->• /. Then, V ;c • (fV x log /) = A x f and d g {fd g f) = d 2 e f 
and we obtain: 

i(V x U)f) - A d e {{d e U)f) - i4V x ■ F x - A id e F 2 - dfiAJ - d\d 2 f = 0, (A.2) 
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where, 


r jo 

Fi(x,9) = / ( gV x V)(x 1 ,9 l ,£ 1 ,x 2 ,9 2 ,£ 2 )dx 2 — -d£ 1 d£ 2 , 


7 r 


F 2 (x,9) = / ((gde 1 V)(x 1 ,9 1 ,£ 1 ,x 2 ,9 2 ,£ 2 ) 


d6o 


+ g(xi, 9 1 , £ 1 , x 2 ,9 2 ,£ 2 )d Bl b(9 1 ,9 2 ))dx 2 — d£\d£ 2 . 


IT 


A.2 Evolution equation for the fiber links 


Following the same principle as for f N and given that the links are maintained over time, 


dr. 


dr. 


i.e. 


*(fc) _ J'(fc) _ 


dt 


dt 


— 0, V k e [1, K], one can write: 


d is K , n = ± E r 


dt 


k =1 




k , dXifk) ,T,/r.k \ dXj^k) 


dt 


+ V„W(Cj‘ wm ) ■ 


V ra *(C‘ 


) • ^37^ + V Xa tt(G* 


<(fc ) J '( fc ) / dt 

k \ ^ 'i(k) <Tf ( 


dt 

dX) 


X 2 ■*v^j(k),i(k)r dt 

d9j[k) 


+ <%i 'H c i(k),m)- d f- + ^i^(Ci(fc),i(fc)) dt 

+ 9e 2 'h(C'i( fc )j( fc )) + de 2 


dt 


= E 1 


E d 


(A.3) 


where Ek corresponds to the A;-th line of (1A.31) . For the sake of simplicity, the computation 
of E\ only is developed here. The computation of the other ones are similar and omitted. 
From Eqs. (12.21) . (12.31) . one obtains: 


E\ 


1 


2K 




dt 


E 

K =1 L 
K 

(Cfttwn)) - V, (u + d log f N )(x m ,e m ) 


k =1 L 


+ v x ,v(cf (i)m ). v x (u + diogf N )(x m .e m ) 


K 


k '=1 


"t ' ^ ^ + V X2 Vd(j(k'),i(k)) j i{k'),j{k’)) 


k '=1 


9 W 1 ^ f (C : ,'(A;) l j(fc)) • I VxiVfi(i(k'),j(k)) + W 2 E(5(j(fc/) J (fc)) j ipi(k'),j(k') 
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where we write V = V(xi, @i, £ 1 , x 2 , d 2 , £ 2 )- Now, exchanging the sums in k and k! and 
using the symmetry of V, one obtains: 

E 1 = - f i((g K ,V xl ^(x 1 ,e 1 J 1 ,x 2 ,9 2 ,£ 2 )-(V x U + dlogf)(x 1 ,9 1 ))) 

K / K 

~{k E ■ £( v ».*(d(*)j<*))««*wo) 

k '=1 ' k =1 

K / K 

~{k E ■ D v »*<^w<»>)4w**’>> 

fc'=l ^ fc=l 

Because there is no restriction on the number of links per fiber, the sums over k cannot 
be simplified in this case. In order to express the third and fourth terms, the number Cf 
(resp. Cj ) of fibers linked to fiber i[k') (resp. j(k')) is introduced: 

Cf = Card({fc | i(k) = i(k') or j(k) = i[k ')}, 

Cf = Card ({k | i(k) = j(k') or j(k) = j(k')}, 

where Card denote the cardinal of a set. Then, as K —> oo, the following expression holds 
for any chosen fiber k'\ 



1 

2 Cf 


K 

k=\ 


—y 

K —^OO 



— dh(U 2 , 


where 


P(Xi( k i), 0i( k i), 0 2 , £ 2 ) = 


g(Xi( k '), 0i( k i),£, x 2 ,0 2 , i 2 ) 


f g(X. 6> i(fc /), £ h x 2 , 0 2 , £ 2 )dx 2 ^-d£id£ 2 ’ 


is the conditional probability of finding a link conditioned on the fact that one of the fibers 
of this link has the same location and orientation as i(k'). Then, as N —» 00 , K —> 00 
such that —> £ > 0 , C* is the mean number of links per fiber. The mean number 
of links in the volume dX^^dOi^ is K f g(X^ k i), 0i( k >), £, x 2 , 0 2 , £ 2 )dx 2 —^d£\d£ 2 and the 
mean number of fibers in dX^^dOi^) is Nf(X i r k /\, 0^'))- Thus: 


Cf 


—y 

N—t 00 
K—to o 


/ g(X i{k , ) ,9i {k/) ,£ 1 ,X2 1 92,£ 2 )dx2^d£ 1 d£2 
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So, we get: 


K 


(^( C m,j{k)) s m,i(k') + ^( c m,%(k)) s m,iO^) 

2e 


fc=l 


->■ 


d0 2 


f( X i(k'),0i(k')) 


K —> oo 
§^>0 


(%)(X i(fc /), Qnjk>),£ 1, x 2 , 6*2, £ 2 )dx 2 —d£ 1 d£ 2 . 


Inserting these expressions in Eq. flA.41) . one obtains: 


E x -f -nig, V xi ^{x 1 ,6 1 ,£ 1 ,x 2l 9 2 ,£ 2 ) ■ (V X U+ dlogf)(x 1 ,9i))) 

JV->oo 
K—*r oo 

f^e>o 

£ K ( 

~^ L ~2K xiV(Ci(k'),j(k')) ' ^1 (-^i(fe')> ^i(fe')) 

fe'=l ^ 


where, 

V'i^i^i) = 77 1 ^ \ [ (g'Vx 1 ^)(x 1 ,9i,£i,x 2: 9 2 ,£ 2 )dx 2 —d£ 1 d£ 2 . (A.5) 

j{Xi 1 0 l ) J IT 

Finally, we find: 


Ei ->■ - p((fif, V. Tl 1 h(a:i, 0i, ^i, ar 2 , 02,40 • (V x f/ + dlogf)(x 1 ,9 1 ))) 

7V-loo ' 7 

K—too 

- ^|3,V I1 V r (ii,fii,^i,i2,^2,4) • V’l(ah, 01)». 


After the same treatment for the four other terms of Eq. (IA.3D and in the limit K,N—> 
oo, % —•> £ >0; one obtains the final equation for g (writting X for (aq, 0i, £\, x 2 , d 2 , £ 2 ))\ 


±((g(X)MX))) 

= -MX), V X1 *{X) ■ w x u(x 1 ,0i)» - v xa tt(x) • X x u(x 2 , e 2 ))) 

- X((g, d ei y(X)d e U(x 1 , e 1 ))) - X((g, d e ^(X)d e U(x 2 , d 2 ))) 

- d^(g, V X1 ^{X) • V. log f( Xl , e 1 ))) - dii((g, V* 2 4/(X) • V, log f(x 2 , 0 2 ))) (A-6) 

- dX((g,d ei ^(X)d e \ogf(x 1 ,9 1 ))) - dX((g, d 02 ^(X)d e log/(x 2 , 6 2 ))) 

- Mg, v,,E(x) • 0i)» - Mg, v^E(x) • M, 0i)» 

- ^£{{g, ipQiV{X) + 9 ^ 6 ( 01 , 9 2 ))xi{ x \, 0i))) 

~ ^{{g, ipe,V{X) + do 1 b(9i,6 2 ))x2{xi, 0i))), 
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where, 


02 (^ 1 , 9 \) = 


/O i,9i 



d9 2 


X 2 


'h) (x 2 , do, £ 2 , xi, 9i,£i)dx 2 — d£ x d£ 2 , (A.7) 


7r 


XifoiA) = J 7 1 , 1 Q y f ( gdo 1 ^)(x 1 ,0 1 ,£ 1 ,x 2 ,0 2 ,l 2 )dx 2 < ^d£ 1 d £ 2 , (A.8) 

X2 (^i, 6 > i) = ^ j (gde 2 ty){x 2 ,6 2 ,£ 2 ,x 1 ,6 1 ,£ 1 )dx 2 ^ L d£id£ 2 . (A.9) 

We introduce the notation Y\ = {xi,9\,£f) and Y 2 = (x 2 ,9 2 ,£ 2 ), and prove the following 
lemma: 

Lemma A.l. For any function h(Yi,Y 2 ), we have: 

((g,h(YuY2)Mxi,0i))) = -((V,. (g(X) 9FYE ),<t(X))), 


{{g,h( Y u Y 2 )M*u<h))) = ~({Vt; (g(X) 

{{g,h(Y u Y 2 ) X i(*i,0i)» = -((5 #I (9(X)AhAii),(p(x))), 

{{ 9 ,/i(U,i2)X2(ll,0l))) = ~((5fe (g(V) 


(A.10) 


f(x 2 ,e 2 ) 

where 0i, 0 2 , Xi and \ 2 are defined by Eq. (IA.5jl and Eqs. (1A.7[1 - (1A.9j) . and where : 


f d6 

F h (x 1 ,9i)= / (gh)(x 1 ,9 1 ,£i,x 2 ,9 2l £ 2 )dx 2 — -d£ 2 d £i 


7r 


(A.ll) 


Proof. Note that for any function h(Y\,Y 2 ), we have: 


((9,h( Y i,Y 2 )'tfi(x 1 ,9 1 ))) 

1 



i) 


>2) 


/(a: i,#i) 


* jo 

(gV xl 'H)(x 1 ,9 1 ,£ /i ,x 3 ,9 3 ,£ 3 )dx 3 —d£ 4 d£ 3 

7T 

dx\ — -d£idx 2 —-d£ 2 

7T 7T 


/ ( /(-r - 1 0 y J ( 9 h)(Y h Y 2 )dx 2 ^-d£ 2 d£ij (g J V Xl ^)(x l ,9 1 ,£ 4 ,x 3 ,9 3 ,£ 3 ) 


dx i — -d£±dx 3 —-df?3 

7T 7T 


= - / V ai ( g(0i,fli,l 4 ,x 3 ,fl 3 ,£ 3 )£ 4 , 373 ,^ 3 , 4 ) 


f{xi,9i 


= -« v “ (9m TO7 )'* m>> 


dx\ — -d£ 4: dx 3 —-df? 3 

7T 7T 
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with Fh defined by (1A.11D . Similarly, we have: 


« g,h{Y 1 ,Y 2 )ip 2 (x 1 ,d 1 )» 
,^2)' 1 





- dQ 

(gV X 2 ^)(x 3 , 0 3 ,4, aq, 0i, £ A )dx 3 —di A di 3 

TT 


doq— -d£idx 2 —-df 2 

7T 7T 


f(x l,0l) 


-^-y J (gh)(Y 1 ,Y 2 )dx2^ L d£ 2 d£^j(g\/ X2 ^)(x 3 ,d 3 ,£ 3 ,x 1 ,9 1 ,£4) 



dx 1 — -d£±dx 3 —^-d £ 3 

TT TT 

dd' 


7)77nOTTr-STT / ^4,4,03, £' 3 )dx' 3 ^d£' 3 d £' 4 

J \ x 2i ^2! J 


7r 


dx\ —- d£[ dx 2 —- d £' 2 

TT TT 


= -J v 4 (gtn, n') j ^§ ) *«. Yi)dx[^-M[dx' 2 ^de ' 2 
= - {{v ’Y (x) lt§£* (x))) 

After the same computations for yq and y; 2 , we obtain Eqs. flA.IOp . 


□ 


Now, lemma [Adi allows us to write the formal limit K, N —>■ 00, 7 —> £ of Eq. (IA.6D 
which reads: 

• (gV a i/(:Ei,0i)+g ^ 0 i) F 1 (o; 1 ,0 1 )) 


— ^de 1 (gdgU(xi, 0i) + £ 


f(x 1,01 


77*1,00) 


X2 ' (9^ xU (x 2 , 02 ) + £ 

-A<% 2 (#d e [/(>2,02)+£ 


f(x 2 , 02 ) 


(072,02)) 


/(ar 2 ,02 


-^2^2,02)) 


(A.12) 


d/i V Xl • ( 


/7i7i 


-VJ(n,0i)) - d/iV 


3,2 


7(^2, 0 2 ) 


v x /(x 2 ,0 2 )) 


-d\d dl ( 9 - dgf(xi,9i)) - dA(9g 2 ( ^ - . dof{x 2 ,9 2 )) = 0, 


77i,di 


7^2,02 


where Ef and E 2 read: 


d0o 


Fi(xi,0i) = / W xl V(x 1 ,9 1 ,£ 1 ,x 2 ,92,£2)g(xi,9 1 ,£ 1 ,X2,92,£2)dx2—d£ 1 d£ 2 , 


F 2 (x 1 ,9 1 ) = 



yV + d 01 b ) ) (xi,9 1 ,£ 1 ,X2,9 2 ,£ 2 )dx2^-d£ 1 d£2- 


Finally, the link creation/deletion Poisson processes, of frequencies 12/ and rq respec¬ 
tively, classically lead to a source term S^g) for Eq. (1A.12D . We recall that a link between 
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two fibers is formed only if the fibers intersect each other, whereas the link deletion process 
obviously acts on existing links only. This leads to the following source term: 

S(g)(x i, 6 >i, x 2 , 6 * 2 , l-i) =Vff{x i, 9 1 )f(x 2 ,9 2 )S(£ i, £(aq, 6>i, x 2 , 6> 2 )6(^ 2 , £(x 2 ,0 2 , aq,6>i) 

- v d g(x 1 ,9 1 ,£ 1 ,x 2 ,9 2 ,£ 2 ), 


where the first term corresponds to the link creation process while the second one, 
to the link deletion process. Here, the quantity f(xi,9i)f(x 2 ,9 2 )S(£i,£(xi,9i,x 2 ,9 2 ) 
S(£ 2 , £(x 2 , 6*2, xi, 9\)dxi^-dx 2 —dl\dl 2 gives the probability of finding a fiber located within 
a volume dx about (aq, 9\) and a fiber located within a volume dx 2 — about (aq, 9 2 ), 
such that they intersect with associated lengths within a volume d£\d £ 2 about (£i,£ 2 ). 
The link creation process generates a new link distribution function proportional to this 
probability at a rate Uf. The quantity —v d g(xi,9i,£i,x 2 ,9 2 ,£ 2 ) corresponds to the decay 
of the link distribution function with rate v d due to the link deletion process. 


B Computation of the non linear term f dg(G[pM]pM)tyd0 

This section is devoted to the computation of the term X 3 given by (15.28(1 . For the sake 
of clarity, the following notations are introduced: 

M = Mg 0 , s 0 = sin 2(9 — 9 0 ), Co = cos2 (9 — 9 0 ). (B.l) 

By symmetry, (h(2(9 — 9 0 ))) = 0 for all odd functions h on [— |, |], where (•) is the average 
defined in Theorem 15.71 We also note from Eq. (14.9p . Hypothesis 16.II and Proposition 15.61 
that we have: 

ArdL 2 rdL 2 


Co = 


aL A 7 


A 8 r)f £48 pc(r) 12£pc(r 

Using Green’s formula, Eqs. (14.8ft , (15.33ft and the same arguments as for X 2 , we get: 


(B.2) 


X 3 = - ( G[pM}pM)d e A> 


d9 


7T 


= -Co 


d9'\ 

7 2 x {pM{ff)) : B(9,9') — )pM(9)d e V 


7r / 


d9 

71 


" 2 2 
7T 
2 


- P c,j \7l(pM(6'))'. (j B(9,ff)M(9)dX^y^- 


P X J XUpM(O')) : J B(9,9')(M(9) - 


Let us first compute V^(pM). We have: 

V x (pM) = MV x p + V X M V x p + V x p ® V X M + pV x M , 
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where V X M is given by (15.301) . A direct computation gives: 


VzM = 2 rM 


soV^o + 2(rso - c 0 )W x 9 0 0 V x 6> 0 


and thus: 


Vl(pM) =M 


v Ip + 2prs 0 V% 


2 rs 0 (V x 9 0 0 V x p + V x p 0 V x 6> 0 ) + 4pr(rsg - c 0 )V x 9 0 0 V x 6> 0 


We now turn towards the computation of 


7T 


2 



7T 


2 


1 d0 


where B(9,9') is given by (I4.10I) . For this purpose, we decompose: 

uj = (uj.ujo)ujo + (tu.cu^a;^ = cos(0 — 9 0 )co 0 + sin(# — 9 0 )uj^, 

where ojq = uj(9 0 ) and such that (u^ay}) is a direct ortho-normal basis of M 2 . Using 
basic trigonometric formulae, one notes that: 


UJ 0 UJ 


1 

2 

1 

2 


(1 + Co)(cUo 0 CUo) + (1 — Co)cj^“ 0 UJq + So[cUo 0 0 CUo] 

I + C 0 [uj 0 0 UJ 0 — UJq 0 CJ^"] + So[cJo 0 UJ o’ + 0 (Uq] , 


where / is the identity matrix. Denoting c 0 = c 0 (9), s 0 = s o (0), Cg = c 0 (9') and s' 0 = Sq(9'), 
we get: 


B(9, 9') = sin 2(9 — 9')[uj 0 uj + uj' 0 a;'] 

1 


— o [ S ° C 0 S 0 C o] 


2/ + (cq + Cq) [cUq 0 UJq — UJ o 0 CJq ] 


(so + Sq) [w 0 0 + Uq 0 cjq] 


= [s 0 c' 0 - s'qCq]I 

+ ^[coSoCq + s 0 Co(6>') - s' o cg(0) - s o c o c o] [w 0 0 UJ 0 - 0 Wq ] 

+ 2 [ s 0 c 0 + S 0 S 0 C 0 — S 0 C 0 S 0 — So Co] [ujQ 0 UJjj~ + UJq 0 (Uo]. 
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Note that B is anti-symmetric, i.e. B(6',6) = —B(9,6'). From the properties of M, we 
get: 

7T 

f d 6 s' d 

/ B(Q,Q ')— = -[w o 0w o ~ u o ®^o‘]-r + + <4®wo]4, 

7 r 4 4 


df) 

M(0)B(6,6')— = —s' 0 (c 0 )I 

71 


^[cUo<8)CUo CUq ® cu 0 ](c 0 Sq(co) + Sq(c 0 )) 
+ \ [w 0 <8) <4 + 4 0 Wo] (Co (So) - So (c 0 )). 


Then, we have: 



1 \d9 

= ~ s ° {Co)I 


+ [wq 0 U1 o — 4 0 44 + [wq 0 4" + W^“ 0 Wo]T2, 


with 


Ti = 


s'o c's'(co) + s(,(c§) 


7 . _ c'o( s o) - 4 (co) c'o 

— — 


4Z 2 2 2 40 2 

Note that this expression is decomposed into an even function T 2 of 6 ' and an odd function 
of 6 ' composed of Sq(co) and 7\. Therefore, (h,Ti) = 0 for all even functions h and 
( h,T 2 ) = 0 for all odd functions h. Moreover, from integration by parts, the following 
relations hold: 


<4 = 


(Co) 


<4 = 1 


(Co) 


( c o) — (Co) - - + 2 


(c 0 ) 


( c o) = l — 2 


( c o) 


4-6 


(B.3) 


(Co) 


M ) = ;(1 - 2 ^ 1 ), 
( 4 ) = 4 ( 1 - 2 % 1 ). 
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Then, 




1 „d0 \dff 


Z 2 Tl ) 7T 


= V^o : [w 0 <8> Uq + cud ® w 0 ] (^2) 

+ 4prV x 6> 0 ® V x 6> 0 : [w 0 ® wjj- + wj[ ® w 0 ] { r(slT 2 ) - (c 0 T 2 ) 


+ 2 r(\/ x p <8> V x 6 0 + V X 0 O ® V x p) : 


(c 0 ) (Sq)/ + [wo ® w 0 - w^ ® w^] (s 0 Ti) 


+ 2pr'V x (y x 9 0 ) : 


- (so)(c 0 )/ + [w 0 <g> w 0 - Wq <g> w^](s 0 Ti) 


(B.4) 


where (using Eqs. (IB.31) and integration by parts): 

(T 2 ) = (rn) 


AZ 2 ’ 

/„ rp \ (c 0 ) 2 1 n (Co). 

( Co 2 } 2r 2 4Z2 ( r )> 

(soT 2 ) = (-- + 2 zjr)jz2 + _ 

(s oTl ) = -M + 3^! + 


2r 2 4rZ 2 ’ 

/ 2rp \ / rp \ ( C o) r 1 1 , 3(Co)-| 

r(s 0 T 2 ) - (c 0 T 2 ) = 1 + — J. 

Then, after some computations and using Eq. (IB.21) . Eq. (IB.41) simplifies into: 


X, = 


dL 2 


24fc(r) 


2 

V 2 (pM(0')) : ( |(B(6»,0(M(6I) - 


= “2^(0 ( " V ' P : [W ° ® Wo± + Wo± 0 ^A^ 

+ Ap(c 0 )V x 6 0 ® V x 6* 0 : [w 0 ® Wq + w^ ® w 0 ](-4y - 1 + ^ 77 ^) 


'AZ 2 

+ 2(c 0 )(pV x V x 6 l o + V x p ® V x 0o + V x $o ® V x p) : 
+ [w 0 ® w 0 - w^- ® w^-] (^ - 1 + ^7^) 

We note that (c 0 ) = c(r). Eq. (IB.5[) leads to (15.351) . 


2 r 
(co)/ 


(B.5) 
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